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Preface 


The  reports  written  by  L.C.  Woods  for  the  Institute  for  Advanced  Technology  (I AT), 
University  of  Texas  at  Austin,  and  collected  together  here  are: 

1.  “The  Current  Melt-Wave  Model,”  October,  1995 

2.  “The  Contact-Spot  Model  of  Transition,”  December,  1995 

3.  “Paper  for  the  8th  EML  Symposium,”  January,  1996 

4.  “An  Armature/Rail  Instability,”  February,  1996 

5.  “Three-Dimensional  Effects  on  the  Melt-Wave  Model,”  May,  1996 

6.  “Boundary  Conditions  Across  a  Rail/Armature  Contact  Surface,”  June,  1996 

A  brief  survey  of  their  contents  and  relationships  follow. 

1.  The  research  started  with  a  rather  idealised  model  of  the  interaction  between  a  rail  and  a 
solid  armature.  It  is  a  two-dimensional  model  and  the  electrical  contact  between  the 
surfaces  is  assumed  to  be  ideal.  The  main  result  is  an  expression  for  the  velocity  vm  at 
which  a  melt  wave  advances  along  the  armature  interface.  It  was  based  on  the  assumption 
that  the  depth  h,  of  the  molten  layer,  assumed  to  be  removed  by  the  viscous  force  due  to 
the  rail,  is  equal  to  the  scale-length  of  the  current  density.  Alternatively,  one  might  have 
taken  h  to  be  equal  to  the  thermal  skin  depth.  Short  of  a  complete  numerical  calculation,  it 
was  difficult  to  decide  between  these  two  possibilities.  Transition  was  assumed  to  occur 
when  the  melt  wave  reached  the  end  of  the  armature  surface  in  contact  with  the  rail.  The 
merit  of  an  ideal  model,  besides  its  simplicity,  is  that  it  gives  a  datum  against  which 
observations  and  numerical  modeling  can  be  compared. 

2.  The  second  report  removes  some  of  the  assumptions  adopted  in  (1),  and  also  replaces 
the  ideal  model  of  electrical  contact  by  the  contact-spot  model.  The  latter  assumes  that  the 
electricity  flows  across  the  interface  between  the  rail  and  armature  through  small,  isolated 
spots,  known  as  “a-spots.”  The  adiabatic  assumption  made  by  earlier  workers  was 
replaced  by  a  diathermic  model  in  which  heat  transport  played  an  important  role.  Transition 
was  deemed  to  have  occurred  when  all  the  spots  were  “burnt  out.”  One  important  variable 
that  was  difficult  to  determine  is  e,  the  ratio  of  the  conducting  area  of  the  spots  to  the  actual 
area  of  the  armature  face. 

3.  This  report  was  largely  a  revised  version  of  (1),  prepared  for  the  8th  EML  Symposium 
in  Baltimore,  held  last  April.  The  main  change  was  that  the  dimensional  arguments  used  in 
(1)  to  obtain  the  scale  lengths  for  the  current  distributions,  were  replaced  by  a  mathematical 
solution.  This  was  obtained  by  using  rectangular  geometry  for  the  armature,  which  is  not 
a  good  representation  of  the  actual  shapes  employed.  However,  as  the  main  purpose  was 
simply  to  provide  a  more  convincing  derivation  of  the  scale-lengths  (valid  in  any 
geometry),  this  was  not  a  real  restriction.  Transition  was  defined  to  occur  when  the  melt 
wave  had  removed  sufficient  material  to  eliminate  the  normal  force  between  the  rail  and  the 
armature.  This  was  more  reasonable  than  the  definition  adopted  in  (1)  and  also  had  the 
advantage  that  the  uncertain  depth  h  no  longer  played  a  role,  since  transition  now  depended 
on  the  product  h  which  is  independent  of  the  mechanism  assumed  to  be  responsible  for 
hv,„  alone. 
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4.  The  movement  of  a  load  along  a  beam,  supported  on  an  elastic  foundation,  can  excite  the 
beam  into  transverse  vibrations  that  have  a  large  amplitude  when  the  speed  of  the  load  is 
close  to  the  velocity  of  free  waves  along  the  beam.  In  this  report,  This  idea  was  applied  to 
the  rail/armature  combination,  with  a  view  to  finding  an  explanation  of  the  critical  velocity 
at  which  rail  gouging  was  likely  to  begin.  The  main  problem  was  to  determine  the 
response  of  the  elastic  foundation  to  a  distributed  load.  For  aluminum  armatures  a  critical 
velocity  of  about  1 ,400  m/s  was  found,  which  is  close  to  that  observed.  This  occurred 
with  symmetrical  displacements  of  the  rails.  With  the  much  heavier  copper  armatures,  to 
obtain  results  close  to  experimental  values,  it  was  necessary  to  assume  the  displacements  to 
be  axisymmetrical.  This  gave  value  of  the  critical  velocity  of  about  600  m/s,  also  close  to 
experiment.  This  work  was  aimed  at  explaining  the  onset  of  gouging,  not  the  details  of  the 
actual  process. 

5.  Because  the  electric  current  has  a  tendency  to  flow  over  the  edges  of  the  armature,  the 
two-dimensional  assumption  fails  at  these  edges.  The  result  is  that  the  melt  wave  travels 
much  faster  near  these  boundaries.  The  report  contains  an  attempt  to  obtain  a  simple 
physical  model  of  the  phenomenon,  avoiding  the  expense  of  an  extensive  computer 
calculation.  The  main  result  is  an  expression  for  the  speed  of  the  melt  wave  as  a  function 
of  the  distance  measured  across  the  armature.  The  report  concluded  with  a  description  of 
how  one  might  set  up  an  efficient  numerical  treatment  that  allows  for  both  three- 
dimensional  effects  and  time  dependence. 

6.  In  this  report  it  is  argued  that  the  current  is  carried  though  the  interface,  mainly  by  the 
penetration  of  the  copper  asperities  into  the  softer  aluminum.  The  relative  motion  of  the 
surfaces  allows  the  asperities  to  “plough”  the  aluminum,  which  action  generates  conducting 
streaks,  that  merge  into  a  continuous  melt  wave.  The  main  change  to  the  melt  wave  model 
of  (1)  is  that  the  expression  for  the  electrical  resistivity  is  somewhat  larger,  since  it  takes 
into  account  some  features  of  the  spot  model.  The  report  concluded  with  an  explanation  of 
the  transverse  striations  that  have  been  observed  with  aluminum  rails.  They  are  attributed 
to  the  instability  of  the  vortex  sheet  formed  by  the  melted  material  between  the  sliding 
surfaces. 

This  work  was  supported  by  the  U.S.  Army  Research  Laboratory  (ARL)  under  contract 
DA  AA2 1 -93-C-0 101. 


Professor  L.C.  Woods 
August,  1996 
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The  Current  Melt-Wave  Model 


Summary 


ll  is  generally  accepted  t.liat,  the  velocity  skin  effect  concentrates  the  current  cross¬ 
ing  the  armature/ rail  interface  at  the  rear  corners  of  a  solid  armature.  Oliniic  heating 
melts  these  corners  and  since  the  resisitivity  of  the  liquid  phase  exceeds  that  of  the  solid 
phase,  the  current  is  deflected  forwards  into  the  solid  region.  The  mechanism  results  in 
a  ‘melt-wave*  moving  from  the  hack  to  the  front  of  the  armature.  There  are  two  distinct, 
descriptions  lor  this  process,  the  entrainment  model  and  the  liquid  model.  In  the  first  of 
these,  shown  in  Fig.  1,  the  molten  metal  is  quickly  removed  by  entrainment  on  the  rails. 
This  leaves  a  high  resistance  vapour  or  plasma  gap  between  the  armature  and  the  rails.  In 
the  liquid  model,  shown  in  Fig.  2.  the  armature/rail  interface  is  initially  filled  with  molten 
metal,  with  vaporisation  occurring  at  a  later  stage.  Since  the  completion  of  vaporisation 
along  the  interface  greatly  increases  the  resistivity,  the  efficiency  of  the  rail-gun  depends 
on  delaying  this  process  as  much  as  possible.  It  is  therefore  important  to  choose  a  model 
of  the  melt-wave  that  accurately  represents  the  physical  situation.  The  speed  vm  at  which 
the  melt  wave  advances  relative  to  the  armature,  is  a  key  parameter.  If  vrn  is  relatively 
large,  then  the  transition  from  metal- to-metal  contact  to  a  plasma  contact  layer  occurs 
early  and  the  velocity  V  attainable  by  the  armature  is  restricted. 


An  expression  for  vm  can  be  deduced  from  the  energy  equation  applied  to  the  solid 
region  of  the  armature  extending  in  front  of  the  melt  wave.  We  shall  show  that  it  is  related 
to  V  by 


Va  IMJ2V 
4 Ijr  pQw2 


(1) 


WheiC  "K  =  ^f§h,{Tm  “  T°)l/’  Q  =  ~  To)  +  Lrn  .  (2) 

In  these  expressions  rja  is  an  average  resistivity  in  the  solid  armature,  r\r  is  the  rail  re¬ 
sistivity,  //,{}  is  the  permeability  of  free  space,  /  is  the  total  current  flowing  through  the 
armature,  p  is  the  armature  density,  w  is  the  depth  of  the  surface  of  the  armature  in 
contact  with  a  rail,  Tm  is  the  melting  temperature,  T0  is  the  laboratory  temperature,  Lrn 
is  the  enthalpy  of  fusion  (latent  heat),  cv  is  the  specific  heat  at  constant  pressure  and  Ka 
is  the  thermal  conductivity  in  the  armature.  In  §5  we  shall  use  (1)  to  derive  a  limiting- 
speed  Vc  for  the  armature. 

To  decide  between  the  models,  we  first  need  to  determine  whether  or  not  the  rate  of 
increase  of  the  volume  of  the  metal  at  inciting  is  greater  or  smaller  than  the  rate  at  which 
it  is  removed  by  entrainment.  For  the  liquid  model  to  be  acceptable,  the  first  of  these 
rates  must  exceed  the  second.  If  this  criterion  is  met,  we  then  need  to  determine  whether 
or  not  the  net  force  acting  on  the  liquid  phase  is  directed  towards  the  armature  or  not.  If' 
not,  then  the  solid  armature  can  accelerate  away  from  its  liquid  phase. 

finally  we  present  a  critical  review  of  work  by  Parks^  on  the  entrainment  model,  and 
by  Janies^  on  the  liquid  model.  Their  expressions  for  vm  disagree  with  each  others  and 
also  with  ours.  We  shall  also  consider  a  work  by  Barber  and  Challita^,  who  introduce1  a 
quite  different  approach  to  the  problem  of  transition,  one  that  challenges  the  validity  of 
either  the  entrainment  nr  liquid  models. 
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1.  Introduction 


Referring  l.o  Figs.  I  and  2.  we  note  I  hat.  in  the  entrainment  model  there  is  no  current 
flowing  across  the  gap  below  CD,  whereas  in  the  liquid  model,  provided  the  conductivity 
of  the  liquid  phase  is  not  much  larger  t  han  that,  of  the  solid  phase,  there  will  be  some 
current  flowing.  Let  R  denote  the  ratio  of  the  resistivity  of  the  liquid  phase  to  the  solid 
phase,  then  R  will  play  a  role  in  determining  the  shape  of  the  region  of  the  liquid  phase, 
shown  in  more  detail  in  Fig.  3.  With  Aluminium  alloy.  R  =  2.4.  whereas  with  Titanium. 
R  =  20.  Therefore  in  the  latter  case  then'  will  be  very  little  current  flowing  across  the 
liquid  and  the  melted  region  will  have  a  shape  more  like  that  shown  in  Fig.  1  for  the 
entrainment  model. 


Figure  1:  Entrainment  model 

First  consider  the  case  illustrated  in  Fig.  2  and  in  more  detail  in  Fig.  3.  As  the  surface 
CD  moves  forward,  increasing  the  volume  of  the  liquid  phase,  either  the  liquid  remains 
with  the  armature,  lubricating  its  motion  along  the  rail,  or  it  is  removed  so  fast  by  the 
friction  force  acting  on  it,  that  a  plasma  gap  is  formed  and  the  liquid  model  is  replaced  by 


Figure  2:  Liquid  model 


the  emtrainment  model.  We!  shall  adopt,  a.  coordinate'  system)  F  fixed  in  the  solid  armature. 
The  rail,  lying  below  the  surface  H,XGX.  move's  with  a  velocity  V  to  the  right.  It  also  has 
ail  acceleration  V  in  t  he  same*  eliremtiem.  Vieweal  from  F.  there'  are  no  ned.  forcers  acting  on 
the  solid  armature,  but  the  liquid  pha.se'  will  be'  subject  to  a  surface  friction  force'  along 
CG.x  and  a  volume  distributexl  electromagnetic  fore:e.  with  a  component;  in  tlm  opposite 
direction.  (Granted  that,  the  me  It- wave  creates  e'nongh  liejuid  to  offset;  the  frictional  losses, 
to  decide  between  the  two  moele'ls.  we  need  estimate's  of  the  magnitude'  of  these  forcers. 

A  related  problem  is  the  distribution  of  the  elee:t,ric  current  in  the  solid  and  liquid 
phases  of  the  armature,  for  this  will  determine'  both  the  rate!  at  which  the  armature  is 
transformed  into  the  liquid  phase!  by  ohmic  heating  and  the  magnitude  of  the!  EM  force! 
tending  to  maintain  the  molten  metal  in  position,  adhering  to  the  solid  armature. 

The  location  ol  the  field  depends  on  the  t  ime  available  for  it  to  diffuse  towards  its 
equilibrium  state.  The  time  7>  for  the  field  to  diffuse  through  a  distance  i  is 

T(:  =  i  =  n/lH)  -  (/'■ o  =  47 r  x  10-' ).  (3) 

where  rj  is  the  resistivity  and  £  is  the  magnetic  (or  resistive)  diffusivity.  And  as  the 
diffusion  time  is  smaller  in  the  melted  material  than  in  the  solid,  once  the  melting  begins, 
the  field  will  diffuse  more  rapidly  through  it.  As  this  happens,  the  EM  force  acting  on 
the  liquid  phase  designated  by  C  in  the  following  -  is  rapidly  reduced.  But  even  if  the 
liquid  phase  had  the  same  resistivity  as  the  solid  phase,  the  penetration  of  the  magnetic 
field  into  the  armature  would  mean  that  the  total  force  acting  on  £  would  be  reduced, 
since  the  elastic  stresses  are  absent  in  C.  For  liquid  Aluminium  alloy  t?  «  16 f2.  i.e.  for 
an  armature  depth  of  10  mm  this  time  is  1.6  ms  and  the  half-depth  is  achieved  in  0.4  ms. 
These  times  are  comparable  with  typical  pulse  times,  so  we  cannot  assume  that;  the  field 
is  in  local  equilibrium.  In  the  solid  phase,  the  diffusion  time  is  about  twice  as  long. 


Figure'  3:  Armature!- Rail  interaction 


2.  Forces  acting  on  the  liquid  phase 

It  is  perhaps  easier  to  understand  the  forces  acting  on  the  molten  metal  if  we  adopt  a 
frame'  /'  that  has  the'  same  velocity  as  F  but  no  acceleration.  In  /•’  the  solid  armature 
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is  accelerating  under  t-lio  EM  Inters  supplied  by  the  current.  II  the  liquid  phase  C  is  to 
’keep  up"  wit  h  t  he  solid  armature.  the  (drees  acting  on  it  must  be  sufficient  t,o  give  it,  the 
same  acceleration.  Hut  we  ea.nnot.  assume  t  ha.t  the  component  of  t  he  magnetic  pressure 
acting  along  the  direction  of  travel  has  the  same  value  on  the  free  surface'  DF ^  a.s  on 
t  Ik'  solid  boundary  DF.  It.  is  more'  helpful  to  consider  the  force'  as  being  t.lie  usual  volume' 
distributed  j  x  B  t.e'rm  that  arts  on  a.  conductor,  with  the'  magnetic  fiedel  B  being  due'  to 
the'  current,  j. 

Consider  the  e:ase  shown  in  Fig.  3.  Suppose  we'  adopt  a  two-dimensional  approximation, 
with  the  magnetic  field  lying  along  OZ  and  the  rail  moving  along  the  OX-axis  in  the 
positive  direction.  In  this  case' 


0D: 

OB- 

(4) 

Ox 

fHYJy,  q  =  (Hl'Jx: 

(j  x  B),. 

(5) 

and 

(j  x  B), 

,  =  -j,Bz  =  -^-(2-B:). 

(0) 

Oy  V  2/j,0  ' ' 

These  equations  give  the  EM  forces  acting  on  C.  There  is  one  other  force  to  consider, 
namely  that  due  to  the  gas  pressure.  But  this  is  such  a  very  small  fraction  of  the  electro¬ 
magnetic  pressure.  B^/2fio  that  wo  can  neglected  it. 

Evidently,  because  of  friction,  some  of  the  liquid  will  adhere  to  the  rail.  As  shown  in 
Fig.  3,  the  thickness  of  this  layer,  b  say,  determines  the  rate  at  which  the  melted  armature 
material  is  removed  from  the  ‘pool’  MCDM.  where  DM  is  a  section  at  the  original  rear 
boundary  of  the  armature.  Let  x  denote  the  distance  from  the  front  edge  of  the  liquid  at  C 
to  M,  then  on  the  assumption  that  the  boundary  layer  on  CM  is  laminar,  its  displacement 
thickness  a  short  distance  downstream  of  DM  is  given  by  the  classical  expression^4) 


where  v  is  the  kinematic  viscosity  and  v  is  the  velocity  of  the  free  surface  DFoo  relative 
to  the  rail.  This  equation  is  based  on  the  assumption  that  the  depth  above  the  layer  is 
infinite,  e.g.  as  in  the  flow  past  a  Hat  plate.  Later  we  shall  verify  that  the  thickness  of  the 
liquid  pool  is  much  greater  than  b\,  which  allows  us  to  set  v  =  V  (see  Fig.  3).  It  follows 
that  the  mass  flux  per  unit  length  of  armature  (in  the  OZ-direct.ion)  is  given  by 

Mass  Flux  —  f)Vb\  =  l.72f){vxV)^~ .  (8) 

In  the  cast'  shown  in  Fig.  1,  the  current  flowing  through  C  is  so  small,  that  in  the 
fra, me  F  there  is  only  the  d  iction  force  acting  on  the  liquid,  whereas  the  solid  phase  is 
accelerated  away  to  the  loft  by  the  EM  forces  acting  on  it.  The  liquid  phase  is  confined 
to  a.  narrow  layer  along  CM.  which  simply  ‘wipes'  the  armature  on  to  the  rail  as  fast  as 
I  hr  melting  occurs. 
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3.  Armature  skin-depths 


TIk1  equation  determining  the  (’volution  of  th(’  magnetic  field  in  an  incompressible  medium 


OB 

—  +  v(.  •  VB  =  B  •  Vv,.  +  £V2B, 

where  v(,  is  l.lie  electron  fluid  velocity, 

Vrx  -j /c.n,,  (rail) 


v,.  = 


-j/c?/.,.  (armature) . 


in  which  x  is  a  unit  vector  parallel  to  the  rails.  Now 

0  _  0 


VB 


r .  o  _  o  i .  _  , 

=  !x*;+y^](,J=^,' 


so  that  from  (4)  we  obtain 

v,  •  VB  = 


-  nojyVz  (rail) 

0  (armature) . 
Thus,  since  B  •  V  =  0,  in  steady  conditions  (9)  yields 

-  /'■  ojyV  (rail) 

0  (armature) . 


(  02  d2. 

+  dyi  1 


(9) 


(10) 


We  shall  use  superscripts  'r  and  V  to  denote  rail  and  armature  values.  By  (4)  we 
can  write  the  equations  as 


dj'x  ,  d3y 


dx 


+  Vfy  =  0. 


Oy 


djl 

dx 


=  0. 


(11) 


It  follows  that  the  length  scale  for  changes  in  j'y  parallel  to  OX  is  8r  =  (,r/V. 

There  are  two  rail/ar mature  boundary  conditions  on  j.  The  first,  follows  from  the 
equation  V  •  j  =  0  integrated  across  the  boundary,  while  the  second  follows  from  the 
x-component  of  Ohm’s  law,  r/j  =  E  +  vxB,  applied  to  each  side  of  the  interface,  and 
the  steady  state  equation  VxE  =  0  integrated  across  the  boundary.  These  boundary 
conditions  are: 


On  y  —  0.  for  all  x  : 


■hi  -hi 


vdx 


Vr3x 


(12) 


Since  /"  and  jy  are  the  same  at  the  interface,  and  is  the  x-length  scale  for  jy,  it 
follows  that  at  least,  just  inside  the  armature,  close  to  the  rail  8r  is  also  the  x- length 
scale  for  y".  Within  the  armature 


Ox  Oy 


=  0. 


Oy 


(13) 


;md  henc<‘  the  scales  lor  spatial  changes  in  each  current  density  component  are  determined 
entirely  by  I, he  boundary  conditions.  The  second  of  ( 13)  shows  that,  j  'j  has  a  //-length  scale 
t  he  same  as  I  lie  c-length  scale  of  /".  i.e.  . 
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I)v  (1 21a 


.&/,  ±(m\ . 

<'1  <)xi)y  <)y  \  Ox  /  ^  <)y  V  Ox  ) 


whore  Ci  =  '/„///•()•  Hence.  different iating  (  1  i)i  by  we 


c 

s« 


0 

Ox 


c 


Ox2 


V 


=  (I. 


on  making  use  of  (12)|  and  (13)_>.  It.  follows  that.  0j"/0y  has  a  T-lengt.h  scale  of  <*>„  =  £„/C. 
Therefore  j“  has  ail  .T-length  scale  of  b„.  From  ( 13)  i  the  .r-length  scale  of;/"  is  the  same 
as  the  y-length  scale  of  j". 

Collecting  these  results,  we  have  for  the  armature  currents  the  following  scales: 


jx  ■  x  1  =  K-  v  1  =  &>■  '■  iw  :  x  T=  <l)r-  y  1  =  > 


with  an  obvious  notation.  We  now  identify  the  thickness  of  the  melt  wave  with  the  y- 
length  scale  of  the  current  density  normal  to  the  amat.ure/rail  interface.  (cf.  Fig.  4). 
Thus 

h  =Sa=UV-  (14) 

The  current  density  normal  to  the  rail/ar mature  interface  can  now  be  expressed  in  the 
form 

:iy  =  n<?lkr(~yli"  (io  =  constant) .  (15) 

The  total  current  passing  through  the  surface  y  —  0  is  therefore 

I  =  w  I  jy  dx  =  'to f)r  jo,  (16) 

7—00 

where  w  is  the  depth  of  the  surface  of  the  armature  in  contact  with  the  rail. 


4.  The  Speed  of  the  Melt- Wave 

Let.  p.  Pi  T  denote  the  local  density,  pressure  and  temperature,  cv  the  specific  heat  at 
constant  volume  and  u  =  cvT  the  specific  internal  energy.  Also  let  TT  denote  the  viscous 
pressure  tensor  and  v  the  fluid  velocity,  then  -TT  :  Vv  is  the  viscous  dissipation.  We  also 
need  the  heat  flux  vector,  q  =  -kVT,  the  energy  input  per  unit  volume  due  to  radiation, 
•</>.  and  the  latent  heat  of  melting,  Lm.  We  shall  assume  that  the  melt- wave  front  is  located 
at  as  shown  in  Fig.  4.  With  these  definitions,  the  energy  equation  becomes^ 

pDu  +  f>pDf>~1  -  -TT  :  Vv  -  V  •  q  +  j  •  (E  +  vxB)  -  pv.,nL.,nt>{xin)  +  •</>,  (17) 

whore  d(xm)  is  the  delta  function,  and 

D  =  +  v  •  V 

<)/. 

is  the  convoctivu  time  derivative. 

Fur  tin*,  solid  pha.se.  we  replace  the  left-hand  sidt'  by  (>\)h  \)j)  ~  pi)//.,  \v i h  is  the 

enthalpy.  C/,71.  and  omit  1.1  H'  viscous  dissipation  and  ra.diaj,ion  terms  from  the  right-hand 
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solid  armature 


Figure  4:  Resistive  skin-depth  in  the  armature 


side.  Then  by  Ohm’s  law, 

pepDT  =  kV2T  +  qj1  -  pLmvmS(x,n).  (18) 

In  order  to  simplify  the  theory,  we  now  assume  that  in  a  frame  fixed  in  the  wave  front, 
conditions  are  steady,  which  allows  us  to  write  DT  =  vmdT/dx  (see  Fig.  4).  Equation 
(18)  now  takes  the  form 

pCpVrn%  =  K  (U? +  ly) +  e2x/dr e~2v/Sa  ~  pL^'nm  ■  (i9) 


We  now  integrate  this  equation  over  the  length  of  the  solid  armature  extending  in 
front  of  the  melt-wave,  which  we  shall  assume  to  be  a  relatively  large  distance.  The 
second  derivatives  of  the  temperature  vanish  at  x  —  —  oo  and  the  temperature  gradient 
in  the  .'^-direction  at  the  melt-wave  front  is  zero.  Thus,  when  (19)  is  integrated  over 
—  oo  <  x  <  0,  the  term  containing  K.d?T/dy 2  vanishes.  Next  we  integrate  over  0  <  y  <  oo. 
In  the  solid  armature,  close  to  the  rail,  we  assume  that  T(x,y)  has  the  form 

T  =  T„  +  (T,n  -  To  )ex'xer!'lx , 


whore  A  is  the  thermal  skin  depth.  On  using  (10)  and  setting  v,n  zero  outside  tin;  range 
0  <  y  <  h.  after  integration  we  arrive  at 


Vn  lk)l2V 

4»/r  (>Qu>2 


-  V,. 


(20) 


W  111 're 


_  /'■()/'■ 

f  n 

(>QVa 


I\))V  -  Q  —  -  71) )  T  Lti 
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5.  The  Maximum  Armature  Speed 


The  acceleration  of  the  armature  projectile  combination  is  given  bv  the  well-known  formula 


(IV  _  vr1 
~dt  ~  ~2M 


(21) 


where  M  is  tlie  total  mass  being  accelerated.  V  is  the  induction  gradient  and  T  is  the 
sum  of  the  various  drag  forces  acting  on  the.  combination.  Eliminating  /  between  (20)  and 
(21).  we  get 


Va  HoM 
4 T/r  pQw2  U 


d.V2 

dl 


+  2VT 


~  Id;  • 


Let.  =  a  bo  the  length  of  the  armature  in  the  OX-direction.  then  integrating  over  the 
time  tm  that  it  takes  the  melt  wave  to  travel  this  distance,  we  get 


a 


V„ 

4  Ijr  pQw2  IJ 


VTdt 


VK  (It  . 


where  Id)  is  the  speed  of  the  armature  at  the  instant  that  the  melt,  wave  is  initiated  and 
Vr  is  its  speed  when  transition  occurs.  Hence 


V2  =  4 rjr  pQvrL1 
Va  Md0 


a  4* 


VTdt. 


(22) 


Notice  that  friction  decreases  Vr.  whereas  heat  loss  to  the  rails,  represented  by  the 
speed  vK.  increases  Vc.  As  we  shall  show  below,  this  heat  loss  is  important:  the  adiabatic 
assumption  adopted  by  some  authors  in  rail-gun  calculations  is  not  justified. 


6.  Application  of  the  Theory 

We  now  return  to  the  problem  discussed  in  the  Introduction.  For  the  liquid  model  to  be 
valid,  we  require  that  the  rate  at  which  the  liquid  metal  is  removed  from  the  pool  MCDM 
(see  Fig.  3)  by  boundary  layer  friction  be  less  than  the  rate  at  which  the  volume  of  the 
pool  is  increased  by  melting.  Thus,  by  (8)  and  (14)  this  condition  is 

(Savm  >  d)  V  =  1.72 [vxVy*  ,  (23) 

where  c  is  the  fractional  volume  increase  due  to  melting. 

Suppose  that  the  armature  metal  is  Aluminium  then  we  need  the  following  physical 
values^ ,  all  expressed  in  SI  units: 

Density:  2700  (2.r>°C),  2040  (000°  C;  solid).  2410  (000  ( h  liquid). 

Enthalpy  of  fusion  3.97x10''.  c;)  =  8.97  x  102.  melting  temperature  T„,  =  000° C. 

Viscosity  1.38  x  10  "'  h  thermal  conductivity  231. 
resistivity  3.7  x  10  (25°  O),  10.7  x  10  ~8  (000UC). 

where  the  resistivity  varies  linearly  with  temperature.  From  these  quantities  we  obtain 
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the  three  dilfusivities: 


thermal  \  (=  Ka/()<),)  =  1.07  x  HP1.  viscous  u  •—  fi.7  x  ]()"“' . 

resistive  6/  (=  Vo /!<■(>)  =  2.91  x  10  '2  (‘JVC).  =  8.51  x  10'-  ((>(»(»“  O). 

With  Copper  rails,  we  also  require  =  1  -ir*  -  10  -  at.  25'  (’. 

For  an  Aluminium  armature  with  Copper  rails,  we  find  from  (7).  (14)  and  (20)  that 

v,„  =  5. 1C  x  Krl(i(//w)2  V  -  9.9  x  K)~'V’  .  (24) 

and  h  =  Sa  =  8.51  x  10“ 2 /V.  #,  =  1.3  x  10_:,(a:/V)5  .  (25) 

To  illustrate  the  theory,  we  shall  apply  it  to  a  rail-gun  experiment  carried  out  in  the 
Institute  for  Advanced  Technology^  'K  Typical  values  are: 

w  =  a  =  25  mm,  /  =  1.2  x  105  A,  V  =  600 m/s.  dV/dt  =  5  x  10r>  m/s2.  In  this  case 
vm.  6.o  m/s,  h  —  142  pm.  t) |  =  8.4  p.m  at  x  =  25  nun.  It  would  take  the  melt-wave  a 
time  tm  =  3.8  ms  to  traverse  the  length  a  of  the  armature  and  initiate  transition. 

More  accurately,  the  relation  between  a  and  t.m  is 
«  =  I<\  jo  I'2V dt  -  I<2  Jtm  V  dt  (K{  =  5.16  x  10~Ui/ve2,  I<2  =  9.9  x  10“  ').  (26) 

Upon  melting,  the  specific  volume  of  Aluminium  increases  by  the  fraction  e.  =  0.088. 
The  melt-wave  therefore  increases  the  volume  of  the  armature  liquid  per  unit  depth  by 
0.088ham  8.1  x  10  in2/s.  This  is  considerably  less  than  the  volume  removed  by  friction, 

which  from  the  left-hand  side  of  (23)  is  5.04  x  10“3  m2/s.  In  this  typical  example,  the 
inequality  in  (23)  fails  to  be  realised  by  two  orders  of  magnitude.  Evidently,  the  liquid 
model  is  not  satisfactory.  The  liquid  near  the  melt-wave  front  will  flow  towards  the  rail 
surface  and  will  be  removed  faster  than  it  can  be  replenished  by  expansion,  with  the 
result  that  a  cavity  will  form  near  the  interface  and  will  extend  along  the  inner  surface 
of  the  liquid  region.  This  is  now  the  entrainment  model  and — according  to  the  above 
theory-  transition  will  occur  as  soon  as  the  melt- wave  reaches  the  front  of  the  armature. 

The  critical  velocity,  i.e.  the  maximum  attainable  before  transition,  is  given  by  (22). 
We  shall  ignore  the  friction  force  T  and  assume  that  the  ‘initial’  velocity  V{)  is  zero.  The 
value  of  L' /M  follows  from  (21).  First  we  shall  omit  the  vH  term.  In  this  case  (22)  yields 
V(:  =  1.45  x  10-1  m/s.  From  V>.  =  tm  dV/dt  we  get  tm  =  2.9  ms,  which  allows  us  to  evaluate 
the  omitted  vK  term.  This  increases  Vr  to  1.57  x  1 0 3  m/s  and  increases  tv,  to  3.1  ms. 

There  is  a.  very  simple  way  of  increasing  Vr.  Armatures  that  are  ‘spring-loaded’  so 
that,  their  surfaces  art*  always  pressed  against  the  rails,  regardless  of  the  loss  of  width 
due  to  ohmic  melting,  would  not  allow  the  vapour  phase  and  the  consequent  increase  in 
the  electrical  resistance,  to  appear.  There  would  be  some  additional  friction  loss  at  the 
beginning  of  the  armature's  motion  relative  to  the  rails,  but  this  might  be  a.  price  worth 
payint*  U>  avoid  transition. 

I3  inally.  tin*  readrr  should  ho  mnindud  that  tin*  throry  of  I. his  roport  and  i.hr  conclusion 
.Pis I'  dra.wn,  is  based  on  a  model  that  may  not  correctly  represent  the  real  physical  sil.na.tion. 
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Experimental  evidence  concerning  the  functional  dependence  ol  Vr  is  required  at  this  stage. 
This  should  be  compared  with  the  dependence  given  by  equation  (22).  Also  the  thickness 
and  uniformity  of  the  Aluminium  deposited  on  the  rails  would  be  a,  useful  guide  to  the 
nature  of  the  melting  process.  To  check  the  point  made  in  the  previous  paragraph,  one 
could  vary  the  initial  rail  pressure  on  the  armature  by  changing  its  size  relative  to  the  gap 
between  the  rails.  If  the  theory  is  correct,  this  should  alter  the  value  of  Vr. 


7.  Some  other  treatments 

Parks*1*  is  concerned  with  developing  a  theory  for  the  entrainment  model.  He  gives  an 
accurate  description  of  the  phenomenon,  but  fails  to  obtain  correct  values  for  the  magnetic 
scale  lengths.  He  assumes  that  the  scale  lengths  are  the  same  in  the  OX  and  OY  directions 
(see  his  Fig.  1).  His  expression  for  h  includes  an  additional  factor  r/ft/r/r  and  consequently 
his  formula  for  vm  is  also  incorrect,  being  3.24(y }r/Jt]a)2  times  the  first  right-hand  term  in 
(20)-  -he  assumes  that  the  process  is  adiabatic.  The  theory  presented  on  the  second  page 
of  his  report  has  little  relevance  to  the  real  problem,  as  the  contradiction  between  his 
figures  1  and  2  immediately  reveals.  He  even  remarks  at  one  stage  that  the  presence  of 
the  gap  between  the  armature  and  the  rails  will  not  appreciably  alter  the  distribution  of 
the  current!  Despite  the  ‘rigorous’  mathematical  analysis,  he  ends  up  by  adopting  rough 
approximations  that  in  effect  reduce  his  approach  to  one  of  dimensional  analysis,  albeit 
wrong.  His  remarkable  agreement  with  experiment  is  fortuitous. 

James*2*  litters  his  work  with  so  many  'form  factors’,  that  it  is  not  easy  to  read,  but  he 
docs  have  a  good  grasp  of  the  physics.  However  he  has  adopted  the  liquid  model  without 
justification.1  Because  he  adopts  the  liquid  model  and  allows  the  current  to  flow  across 
it  (see  Fig.  2),  James  defines  the  current  wave  to  be  a  vaporisation  wave.  To  modify 
equation  (20)  to  apply  to  this  case,  we  need  only  to  increase  the  value  of  Q  to  include  the 
additional  internal  energy  between  the  melting  and  vaporisation  temperatures  plus  the 
heat  of  vaporisation.  Apart  from  some  form  factors,  James’  formula  for  vm  is  similar  to 
that  given  by  the  first  right-hand  term  in  (20). 

Barber  and  Challita*3*  have  no  confidence  in  either  the  liquid  or  entrainment  models, 
which  they  find  irrelevant  because  the  assumption  of  a  perfect  contact  surface  between 
rail  and  armature  is  never  attained  in  practice.  The  electric  current  flows  through  discrete 
spots  that  occupy  only  a  small  fraction  of  the  surface,  and,  they  argue,  this  considerably 
weakens  the  influence  of  the  velocity  skin  effect. 

They  consider  the  problem  discussed  above  about  the  rates  of  production  and  loss  of 
liquid  and  conclude  that 

The  mass  rate  of  material  removal  [due  to  friction]  is  smaller  than  the  mass  rate  of 
expansion  at  all  velocities  of  interest.  The  ratio  is  about  30%  at,  1  krn/s  and  increases 
with  velocities  to  about,  70%  at.  3  km/s.  Therefore  we  do  not  expert,  a  loss  of  lubrication 
and  do  not  expect  solid-solid  friction  to  affect,  contact  transition. 

1  About,  two  years  ago  Trevor  James  asked  me  to  help  him  find  supporting  arguments  for  this  model.  I 
considered  onlv  the  forces  acting  on  the  liquid  and  derided,  with  making  a  serious  calculation,  that  it  was 
plausible.  At  that  time  1  did  not.  consider  the  competition  between  removal  l>v  friction  and  replacement 
by  expansion,  which  1  now  think  settles  I  lie  rase  in  favour  of  tin*  entrainment  model. 
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This  conclusion  appears  to  he  in  conflict  with  the  one  we  deduced  from  (23).  Their 
theory  differs  in  detail,  because  they  ignore  the  latent  heat  of  melting,  hut  the  principle 
is  the  same.  The  two  results  are  reconcilled  when  it  is  realised  that  the  ratio  of  the  right- 
hand  side  of  (23)  to  the  left-hand  side  depends  on  the  dimension  of  the  contact  surface  to 
the  power  1.5.  A  typical  spot  has  a  radius  of  100  //in.  whereas  the  armature  surface  might 
be  10  mm  or  more  across.  It  follows  from  these  figures  that  in  the  spot  model,  the  ratio 
just  defined  is  about  10“,}  smaller  than  in  the  liquid  model.  This  is  sufficient  to  satisfy 
the  inequality  corresponding  to  (23)  and  to  justify  the  conclusion  just  quoted. 

What  remains  to  be  considered  is  the  extent  to  which  the  very  high  relative  motion  of 
the  contact  surfaces  affects  the  spot  model,  which  was  originally  developed  for  stationary 
surfaces.  Experimental  evidence  for  the  model  might  be  a  streaked  appearance  of  the 
Aluminum  deposits  on  the  rails.  We  shall  give  a  detailed  account  of  the  spot  model  in 
another  report. 
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Summary 

In  Progress  Report  no.  I  (hereafter  «I«mk>U*(1  by  [  1  ]),  thorn  is  an  account,  of  the  melt- 
wav(<  model  for  transition  in  solid-armature,  rail-gun  launchers.  This  model  contained 
a  number  of  assumptions  concerning  the  electrical,  mechanical  and  thermal  interactions 
between  the  rail  and  the  armature.  These  are  (i)  perfect  electrical  contact,  (ii)  the  removal 
of  liquid  metal  by  viscosity  and  (iii)  adiabatic  conditions  except  possibly  in  the  region 
upstream  of  the  melt-wave.  As  our  main  objective  was  to  determine  the  speed  of  the  melt- 
wave  relative  to  the  armature  and  to  use  this  to  calculate  the  projectile  speed  at  which 
transition  occurred,  only  the  first  of  these  assumptions  was  critical.  However  to  determine 
the  thickness  of  the  material  deposited  on  the  rail  by  the  passage  of  the  armature,  we 
made  use  of  (ii)  and  (iii).  We  now  believe  that  these  are  not  good  approximations,  at  least, 
for  the  entrainment  version  of  the  melt- wave  model  and  in  §1  we  shall  change  them.  The 
modified  theory  will  then  be  used  to  determine  the  thickness,  Svu  of  the  deposited  layer. 

In  [1]  we  referred  to  the  work  of  Barber  and  ChallitaW  concerning  the  electrical  contact 
between  the  armature  and  the  rail.  Their  model  of  armature  contact  transition,  first 
presented  at  the  5th  EML  Symposium^,  is  based  on  ‘adiabatic  accumulation  of  resistive 
dissipation  in  contact  spots’.  Such  spots  are  usually  known  as  ‘a-spots’1 .  In  the  melt- wave 
model,  as  shown  in  [1],  the  velocity  skin  effect  (VSE)  plays  a  central  role  in  determining 
when  transition  occurs.  However,  with  their  model  of  electrical  contact,  Barber  and 
Ghallita  find  that  VSE  has  little  effect  on  transition.  Their  treatment  has  some  puzzling 
features,  but  their  main  point,  namely  that  the  melt-wave  model  is  poor  representation 
of  the  annature/rail  interaction,  is  convincing.  Our  main  task  in  this  report  is  to  develop 
a  new  theory  for  an  ‘a-spot’,  in  thermal  contact  with  its  environment,  and  to  use  it 
to  determine  the  thickness  of  deposited  material  and  the  time  (of  flight)  at  which  we 
can  expect  transition  to  occur.  The  nature  of  the  deposited  layer  enables  us  to  decide 
between  the  melt-wave  and  a-spot  models,  since  with  the  latter,  we  would  expect  to 
see  a  deposited  layer  of  variable  thickness,  striated  in  the  direction  of  motion.  Some 

preliminary  measurements  of  6m  obtained  at  the  Institute  for  Advanced  Technology!3) 
show  this  feature. 


1.  The  Entrainment  Model  Modified 


In  Figs.  1  and  2  we  show  the  two  versions  of  the  melt- wave  model  that  were  considered  in 
[1].  The  liquid  model  is  valid  only  if  the  rate  at  which  the  volume  of  the  pool  is  increased 
by  melting  exceeds  the  rate  at  which  it  is  removed  by  friction  with  the  rail.  In  this  case 
is  reasonablfi  to  *®uuie  the  deposited  thickness,  6m,  is  approximately  the  same  as 
the  displacement  thickness  of  the  boundary  layer  that  starts  very  dose  to  the  melt- wave 
Iront,  where  the  fluid  is  relatively  deep.  However  the  condition  on  the  rates  is  not  met,  so 
this  model  is  not  satisfactory.  With  the  entrainment  model,  the  fluid  is  never  ‘deep’  so 

the  boundary  layer  is  not  a  relevant  concept.  In  [1]  we  omitted  to  complete  this  model- 
w<*  shall  do  Uial,  hurt;. 

A  more  important  matter  is  the  question  of  the  thickness  of  the  gap  between  the  soli, I 


W«"  111,*  (pi.tini  reason  dial,  in  (lie  old  . . .  ( lie  snl,|ert .  then  radius  w;i 

lei  I  (•«  n 1 


»s  ,itw;i\\s  ildiolril  hy  I  hr 
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armature  and  the  rail  behind  the  advancing  melt-wave,  denoted  by  h  in  Lig.  2.  In  }j3  of 
[1],  still  adhering  to  an  adiabatic,  description  of  tin1  process,  wo  argued  tha,t  this  should 
lx*  equal  to  the  y-Iength  scale*  of  the*  current  density  ma  nia!  to  tin*  armature/rail  interface 
just,  upstream  of  the  melt-wave.  Therefore,  by  the  theory  developed  for  this  scale,  we 
wrote 

h  =  K  (A«  =  6i/^-  s  -  '///'oh  (1) 

where  V  is  the  relative  speed  of  the  rail,  7/  is  the  resistivity  and  (a  the  magnetic  diffusivity 
of  the  armature.  One  defect  of  this  argument  is  that  btl  is  the  y-lcngth  scale  for  the 
current  density  in  the  same  direction,  i.e.  of  jy,  only  close  to  the  interface  boundary  and 
in  front  of  the  melt-wave.  Since  the  components  jx  and  jy  satisfy  Laplace’s  equation  in 
the  armature,  their  scale  lengths  are  determined  by  the  nearest  boundary  condition.  Our 
argument  inverts  this  we  have  allowed  the  scale-length  appropriate  to  one  part  of  the 
boundary  to  determine  the  boundary  condition  for  another  part.  The  scale-length  for 
jy  in  the  direction  parallel  to  the  wave-front,  i.e.  perpendicular  to  the  interface,  is  to  be 
determined  from  h,  not  the  other  way  round.  So  what  does  determine  h? 

The  adiabatic  assumption  is  the  problem,  for  it  must  be  the  thermal  properties  of 
the  armature  that  will  determine  the  range  over  which  the  heat  from  the  intense  electric 
current  will  extend.  Of  course  if  we  were  to  adopt  a  mathematical  approach  and  write 
down  all  the  equations  and  boundary  conditions,  then  no  doubt  the  ‘truth’  of  the  matter 
would  emerge  automatically,  but  this  would  be  a  very  difficult  calculation,  for  behind 
the  melt-wave,  we  have  a  moving  boundary,  whose  location  is  initially  unknown.  Also  we 
would  need  to  include  time  in  the  list  of  independent  variables,  so  increasing  the  dimensions 
to  three— -or  four,  if  the  true  geometry  of  the  armature  is  admitted.  We  shall  avoid  this 
complexity  by  adopting  dimensional  arguments,  but  avoiding  empirical  elements. 

Let  k,  p  and  cv  denote  the  thermal  conductivity,  density  and  specific  heat  at  constant 
pressure  of  the  armature  material,  then  the  thermal  diffusivity  is  defined  by  x  —  ft/pcp. 
This  has  the  dimensions  of  a  length  squared  over  a  time.  Therefore,  if  t  denotes  the  time 
during  which  the  diffusion  takes  place,  the  distance  of  penetration  of  the  heat  is  (*<)*• 
This  distance  determines  the  width  h.  Hence  (1)  is  replaced  by 

h  =  (x<)*  (x  =  K/pCp\  (2) 

Let  a  denote  the  acceleration  of  the  armature,  then 

of  =  ( kL'1 2  -  Fa~  Fd  -  Fp)/(ma  +  mp),  (3) 

where  ma  is  the  armature  mass,  rnp  is  the  projectile  mass,  k  is  a  constant,  equal  to  £  or  a 
little  smaller,  V  is  the  induction  gradient  of  the  gun,  /  is  the  total  gun  current,  Fa  is  the 
ablation  drag,  Fd  is  the  viscous  drag  and  Fp  is  the  retarding  force  due  to  the  background 
gas.  In  a  complete  theory,  we  should  use  (3),  but  to  obtain  approximate  results  for  the 
thickness  of  t  he*  layer  on  the  rails,  wo  shall  bo  content  to  take  <\  from  experiment  and  also 
to  give  it  a  constant  value;  experiment,  indicates  that  this  is  roughly  so  over  its  middle 
range;  of  values. 

The  assumption  of  constant  acceleration  allows  us  to  write 

1  o  / ,  .1  \(\t“.  I  ( 'l  r  / 1  (  j ) 
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Figure  1 :  The  Liquid  Model 


solid  armature 


wher?  V  is  the  armature  velocity  relative  to  the  rails  and  x  is  its  displacement  at  time  t. 
We  can  now  write  (2)  as 

h  =  x*2*(z/a)*  [^4.63x10-4x?],  (5) 

where  the  numerical  expression  in  (5)  is  for  Aluminium,  (x  =  1.07xlO”4),  with  the  (ex¬ 
perimental)  average,  a  =  5><105  ms"2.  (We  are  not  aiming  at  accuracy  at  this  stage,  but 
some  rough  estimates  with  which  we  can  compare  experiment.) 

Now  consider  the  conservation  of  mass  in  a  frame  fixed  in  the  melt- wave*  front.  With 
subscripts  s  and  f,  denoting  the  solid  and  liquid  phases,  and  with  v1n  demoting  the  melt- 
wave  velocity,  on  the  assumption  that  the  mass  entrained  on  the  rail  equals  the  mass 
melted  in  the  wave  front,  we  get 


/V;?nh  -  /'/'O'  1  t}m)^in- 

Hence  the  thickness  of  the  deposited  layer  is  given  by 


A 


m 


I1  !>s 

I  I1  I'! 
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(7) 


liquation  (24 )  of  [  I  ]  y  i<4<  is 

it  ■;  v„,/V  -  5.  lG*(//ni)2  [a  Kx |()  |:’/2], 

wIk'cc  7/j  denotes  i.h('  drpl.h  of  t.lio  arinal.un*  surface  in  contact,  with  (  lie  rail  a, ml  lor  Ui<‘ 
value'  in  brackets.  we  adopted  w  25.4  mm.  As  fi  «  1,  and  (>s/(>(  2(i40/24l9  -=  1.0%. 

equations  (5)  to  (7)  yield 

K,  --=  8.7()xl()-13  I2 h  =  4.05x10" Hi  l2x- i  .  (8) 

For  example,  if  /  =  3.5*i():*  A  and  x  =  50  mm,  then  h  —  219  /im  and  A,n  =  23  /an. 

According  to  this  theory,  the  deposition  thickness  is  a  maximum  where  I2x *  is  a 
maximum,  which  will  occur  a  very  small  distance  beyond  the  station  where  the  current  is 
a  maximum. 

It  is  of  interest  to  compare  this  with  the  boundary  layer  thickness,  namely 

<*>,  =  1.72  (vx/V)*  «[4.11x10"M],  (9) 

using  the  constant  acceleration,  a  =  5x10s  ms-2.  The  thickness  are  similar  for  typical 
currents,  /  «  3><105  A. 

The  characteristic  width  of  the  current  layer  just  in  front  of  the  melt-wave  is 
sr  =  £r/v  =  a {2 ®a)i  «  1.35xlO-r>/.T2  , 

using  the  value  of  a  given  above.  At  x  —  50.8  mm,  this  gives  60/iin,  a  surprisingly  narrow 
strip. 


2.  The  Adiabatic  Spot  Model 

On  the  melt-wave  model,  transition  occurs  when  the  melt-wave  reaches  the  front  of  the 
armature.  Here  we  shall  consider  the  rather  different  description  of  transition  proposed 
by  Barber  and  Challita^2),  based  on  the  adiabatic  accumulation  of  resistive  dissipation 
in  a-spots.  Their  first  step  was  to  use  time-scale  arguments  to  establish  that  there  was 
insufficient  time  for  heat  to  be  lost  from  the  spot  by  thermal  conduction.  We  shall  review 
this  key  phase  of  their  argument. 

We  shall  start  with  a  general  form  for  the  energy  equation  in  a  fluid,  with  a  fluid 
velocity  v  along  OX,  sheared  in  a  direction  along  OY  (see  Fig.  1).  Let  //  denote  the 
kinematic  viscosity  (or  viscous  diffusivity),  then  the  viscous  dissipation  is  2 pu(dvfdy)'2; 
we  also  need  the  resistive  dissipation,  yj2,  due  to  current  density  •/,  the  heat  flux  vector, 
q  —  — /cV7\  and  the  latent  heat  of  melting,  Lm.  In  constant  pressure  changes  the  energy 
equation  reads: 

~  F  VJ1  +  2pi'  -  i>Lmvm6(x,„ ) ,  (1(1) 

() 

I)  v  -  l  V  -  V 


wlietV 
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is  t,h('  cot ivccl.ivi*  l.imc  derivative,  /)(./:,„)  is  the  delU  fund, ion.  and  vw  is  t.i le  velocity  ofl.hr 
melt -front,  al  x  --  xw. 

do  compare  the  Imus,  we  define  time  scales  rv  by  ()cvAT/tx  •-  A',  where  X  denotes 
a  ritflil.-lia.nd  term  in  (10)  and  AT  is  a.  typical  temperature  chantfe.  Let  a  denote  a. 
characteristic  lentfth.  such  as  the  spot  radius,  and  define  the  thermal  diffusivitv  by  y  = 
ti/fH :?M  then 

Tk  =  «7x  ,  T„  =  (>cv AT/ (Ilf  ) ,  r„  =  c]lAT6i/(2,yV2) .  (1 1 ) 

whore  for  the  viscous  length  scale,  we  have  taken  the  boundary  layer  thickness  defined 
in  (9).  To  deal  with  the  last,  term  in  (10).  we  integrate  over  a  typical  length  and  take 
Tm  =  a/vm-  We  adopt  the  ohmic  heating  term  for  comparison.  Thus 

Tm  =  pLm  /  {''ij  )  1  (12) 

gives  the  adiabatic  melting  time  due  to  ohmic  heating.  We  shall  take  AT  to  be  the 
temperature  increment  required  to  raise  the  armature  material  to  its  melting  point,  then 
Tv  +  Tm  is  tll(-  Wtal  time  necessary  to  melt  the  material.  We  shall  find  it  numerically 
convenient  to  replace  j  and  a  by  g  =  10“loy,  and  c  =  104a.  In  SI  units,  we  have 

X  —  1.07*10“4  ,  p  =  2640 ,  AT  =  035,  cv  =  8.97*102  ,  Lm  =  3.97M0r’ , 

r)  =  0.72X10"7,  v  =  5.7x10“' ,  j  =  10105,  a  =  10“1c. 

Using  these  values  and  (9),  we  arrive  at  the  time  scales 

tk  =  0.93X10-V  ,  Trj  =  1.90x10- '/g2  ,  r"  =  lA6^10~A/g2  ,  Tl/  =  8.44X10-1®-1/2  . 

(13) 

Evidently  viscous  heating  is  important  only  at  large  values  of  the  displacement  x.  We 
shall  ignore  it  in  the  following.  With  the  values  of  /,  w  and  6r  given  at  the  end  of  §1, 
we  find  j0  ~  2*10n  A  in-2,  which  suggested  the  factor  1010  in  defining  g.  It  remains  to 
choose  c  and  g. 

The  adiabatic  assumption  is  valid  only  if  tk  »  tv,  i.e.  if  eg  »  1.43.  In  their  first 
paper,  Barber  and  ChallitaO  accept  the  values  c  =  10  and  0.1  <  g  <  1  as  being  typical; 
the  adiabatic  constraint  is  not  satisfied  at  the  lower  currents.  However  in  their  second 
paper*1),  they  reduce  c  to  1  and  have  extended  the  range  for  g  to  20.  It  follows  that  over 
the  current  range  0.1  <  g  <  1.43  ,  spots  of  this  size  are  not  adiabatic.  In  fact,  at  values  of 
a  —  100 /nil  and  j  =  1010  Am  ",  tk  <  r,;,  which  means  that  melting  temperature  cannot 
be  reached.  In  the  steady  state  AT  is  determined  by  tk  =  r,;,  i.e.  AT  =  /k. 

It  is  reasonable  to  assume  that,  as  with  the  entrainment  model  of  Pig.  2,  transition 
occurs  when  the  material  in  physical  contact  at  a  point  has  completely  melted.  Using  the 
adiabatic  hypothesis,  Barber  and  Challita,2)  arrive  at  the  empirical  law. 

fl{  ■> 

k  /  (/"//'’),//  A  (/.:  =  const.),  (14) 

where  /.„  is  the  conduction  starting  time,  is  the  time  of  contact  transition,  A  is  the; 
nominal  contact  area,  k  is  the  normal  force  on  the  contact  and  /,  is  a  const. ant.,  whose 
value  depends  only  on  the  contact  materials.  Thus  transition  can  be  delayed  by  increasing 
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armature 


Figure  3:  A  Contact.  Spot 


the  contact  force,  increasing  the  contact  area,  or  choosing  a  material  with  a  low  transition 
constant  k.  They  found  that  the  ‘constant’  k  was  not  velocity  dependent  up  to  1.2  km/s. 

The  first  two  of  these  conclusions  are  scarcely  surprising  and  apply  equally  to  a  modified 
version  of  the  melt-wave  model.  Evidently  if  there  is  a  normal  force  on  the  contact,  it  will 
not  be  sufficient  for  the  melt-wave  to  reach  the  front  of  the  armature  to  initiate  transition, 
for  at  this  point,  or  perhaps  a  little  earlier,  the  normal  force  will  squeeze  out  the  liquid 
between  the  surfaces  and  re-establish  solid  metal  contact.  This  will  continue  until  the 
normal  force  falls  to  a  low  value. 

In  their  1993  paper,  Barber  and  Challita  extended  their  work  to  greater  speeds.  They 
found  that  the  transition  ‘constant’  remained  almost  unchanged  up  to  speeds  of  2  km/s, 
allowing  them  to  conclude  that  the  velocity  skin  effect  had  little  influence  on  the  results, 
apart  from  shifting  the  current  from  leading  spots  to  trailing  spots. 

3.  The  Diathermic  Spot  Model 

We  shall  now  develop  a  theory  in  which  heat  loses  from  the  a-spot  play  a  significant  role  in 
the  phenomenon  of  contact  currents.  This  is  an  extension  of  the  classical  treatment  of  the 
subject*3).  Fig.  3  illustrates  an  idealised  contact  spot.  It  is  drawn  on  the  assumption  that 
the  armature  is  the  harder  material,  but  this  may  lx;  incorrect.  However  interchanging 
the  rail  and  armature  labels  in  the  figure  will  have  no  effect  on  tlx;  theory  to  follow.  As 
we  have  drawn  it,  initially  an  armature  protuberance  indents  the  rail,  reaching  pressures 
up  to  half  of  that  necessary  to  make  the  indented  material  yield*3).  With  copper,  this  is 
about  100  to  200  MPa..  The  intense  flow  of  electric  current  through  the  restricted  region 
of  contact,  heats  it,  completes  the  softening  process  and  an  a-spot  is  created.  The  central 
region  quickly  reaches  the  melting  temperature  and  with  continued  ohmic  heating,  the 
liquid  region  spreads  outwards.  The  relative  motion  of  the  rail  and  armature  removes  the 
molten  metal  as  fast,  as  it  is  created  and  when  the  spot  is  ‘consumed’  in  this  way,  metallic 
contact  between  the  rail  and  armature  is  lost  at  this  point.  As  a-spots  are  ‘burnt  out’, 
new  spots  will  be  able  to  form,  as  the  surfaces  come  infinitesimally  closed'. 
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Inil.ially  I  lie  <I<'|)(.li  of  indentation.  a  is  determined  by  the  hardness  of  (.lie  copper  j,,  |,|1(. 
rail.  lull,  as  the  heating  and  softening  rapidly  proceeds,  the  thermal  skin  depth  given  by 
(•0)  will  play  a  role.  Which  of  these  processes  is  dominant  remains  to  be  discovered:  for 
om  discussion,  we  shall  adopt  the  thermal  skin  depth  for  the  same  reasons  as  advanced 
in  . 

In  order  to  simplify  the  theory,  we  have  modified  the  geometry  of  the  spot  shown  in 
Fig.  3  to  produce  the  axially  symmetric  representation  shown  in  Fig.  4.  The  figure  contains 
a  central  region  of  molten  metal,  expanding  radially  outwards  with  the  progress  of  time. 
In  drawing  this  figure,  we  have  assumed  that  the  copper  rail  is  harder  than  the  armature 
material.  When  the  melt  wave  reaches  the  spot  boundary  at  r  =  a,  local  transition  will 
occur.  But  with  a  multitude  of  a-spots,  potential  and  actual,  this  should  not  be  confused 
with  armature  transition. 


armature 


Figure  4:  The  Contact  Spot  Model 


We  shall  first  develop  a  theory  for  the  stationary  a-spot.  We  start  from  (10),  omitting 
the  unimportant  viscous  heating  term.  We  integrate  over  the  range  ()</.<  T(,’  where  r0 
is  the  time  at  which  the  a-spot  first  reaches  melting  temperature.  We  cannot  accurately 
represent  the  first  right-hand  term  short  of  a  full  calculation,  with  appropriate  boundary 
conditions.  Ihe  sketch  in  Fig.  3,  shows  that  this  would  be  very  dependent  on  the  details 

of  the  spot  geometry.  We  shall  therefore  be  content  with  the  approximation  adopted  in 
(11).  Let 


(vf)  =  —  f  1  /j'2  (If 
T{ )  ./(I 

then  (10)  yields  r0  1  =  t,~1  -  r“‘,  i.e. 


_  pcpAT 

T’'  ~  (rip)  ' 


(15) 


T() 


T,,Tk 
tk  ~  T,, 


(T<.  >  n,)  • 


Adopting  the  constant  values  of  //  and  j  used  in  (13),  this  gives 


1  .()«)>  10  V 
(<•</)•'  2  01 


('•</  '•  1,13) 


(Hi) 


(17) 
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When  I, he  melting  temperature  has  been  achieved,  the  spot  remains  at  constant  tem¬ 
perature  during  t,h<'  passage  of’ a,  melt  wave  from  the  a.xis  f,o  the  spol.  boundary.  For  this 
case  (10)  yields 


<IT  K  d  (IT 

fx-i^nr—  =  r  f  if;r  ~  f>bmVm(r)b(r). 

dr  r  dr  dr 


We  shall  assume  that  the  cylindrical  contact  region  is  all  at  tin1  melting  temperature,  so 
that  the  left-hand  term  is  zero,  Multiplying  by  lhxr  and  int.egra.ting  over  0  <  r  <  a ,  we 
find 

^7/i  ”  (^  VJ  ~~  ^  ) / [2f)Ij,n 7  ).  (18) 

where  we  have  written  [dT/dr]r~n  as  -  A T/a. 


The  time  required  to  melt  the  whole  of  the  cylinder  is  r,„ 
(11)  and  (12),  we  obtain 

1  1  /2cpAT\  1 

Tin  rw  '  Lm  '  Tk’ 


l()dr/vm.  Hence,  using 


(19) 


Adding  tq  and  rm,  we  arrive  at  the  time  for  spot  transition,  rtsr  : 

tvtk 


T,„  = 
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Tk  -  T?]  rK  ~  2.87 T"  ' 
on  using  the  numerical  values  given  in  §2.  Hence  from  (13), 
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(eg  >  2.12). 


(20) 
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A  typical  value,  with  c  =  1.5  and  g  =  5,  is  rtAr  =  14  /is. 


4.  Modification  Due  to  the  Velocity  Skin  Effect 


It  is  obvious  from  Fig.  3  that  rapid  relative  motion  between  the  armature  and  rail  will 
render  the  above  theory  irrelevant.  For  example  at  V  =  500  ms-1,  the  diameter  of  the 
example  at  the  end  of  the  previous  section,  2a  =  3X10-4  m,  will  be  traversed  in  0.6 /xs, 
long  before  the  melting  temperature  is  achieved. 

The  velocity  skin  effect  will  cause  a  concentration  of  the  electric  current  near  the  rear 
of  a  sliding  contact  spot.  The  heating  will  not  be  symmetric,  as  depicted  in  Fig.  3,  but 
will  result  in  a  melt  wave,  that  passes  through  the  spot  with  a  velocity  vln ,  as  illustrated 
in  Fig.  5.  Equation  (20)  of  [1],  viz. 
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will  apply.  Hence 
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Figure  5:  A  Melt  Wave  Crossing  an  A-Spot 


The  spot  transit  time  is  therefore 


r,  =  0.105 


2aV 
So  ' 


(23) 


The  thickness  6m  of  the  deposited  layer  is  given  by  (6),  modified  for  the  present  appli¬ 
cation.  From  (4)  to  (6), 

ft  =  vm/V  =  9.53(</o/V)2,  h  =  1.23>«10-2(x/a)i,  V  =  (2ax)?  .  (24) 

Hence  6m  =  P(ps/pc)h  =  6.42x10-2-^.  (25) 

There  remains  the  difficult  question  of  defining  transition  for  the  armature  as  a  whole. 
The  velocity  skin  effect  will  concentrate  the  current  into  the  a-spots  closest  to  the  rear  edge 
of  the  armature  face  in  contact  with  the  rail.  To  facilitate  the  description,  imagine  these 
to  be  in  rows  parallel  to  the  trailing  edge.  As  each  row  of  spots  is  eliminated,  the  current 
will  move  forward  to  the  spots  in  the  next  row,  and  so  on  until  the  front  row  is  reached.  Of 
course  it  may  not  happen  in  this  way— new  spots  might  appear  behind  the  rows  already 
visited.  But  let  us  assume  that  row  by  row  elimination  is  a  rough  approximation  to  the 
process. 

There  is  one  important  variable  that  is  particularly  difficult  to  determine.  This  is 
the  ratio  e  of  conducting  area  in  the  spots  to  the  actual  area  of  the  armature  face,  or 
equivalently,  the  number  density  Ns  =  e/(ir a2)  of  spots  of  radius  a.  Values  of  £  of  between 
0.1  and  0.2  have  been  quoted^.  If  we  rearrange  the  circles  into  squares  of  the  same  area, 
we  find  that  the  conducting  length  per  unit  length  of  armature  is  ei.  Then  if  w  is  the 
width  of  the  armature  measured  parallel  to  the  rails,  the  conducting  length  is  we?,  and 
the  transition  time  for  the  armature  is  r„,  where 

r„  =  vie  J  /a,,,  =  ().  105  welv/g%  .  (26) 


To  illustrate  the  theory  with  some  numbers,  let  us  assume  that  a  5  -10''  ms”2 


</o  =  10,  a  —  1 .5*10" 4  m,  w  =  0.025  m,  w  - 
get 

T,.  71  //.s,  F  224  ms”1,  n„,  ~  4.25  ms 


0.05  in,  e  0.1.  Then  from  (23)  to  (2(i)  we 
h  -  218/im,  (i,„  -5 //.in.  r„  -  I. Klims. 
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The  next,  stage  in  the  research  requires  some'  experimental  data  to  allow  us  to  confirm 
or  reject  the  above  theory.  Examination  of  the  surface's  of  the  armature  and  the  rail  after 
rail-gun  firings  should  provide  valuable'  dues.  However,  oven  if  the  theory  is  correct  so  far 
as  the  thermal  effects  arc'  concerned,  transition  may  have'  a  different  cause. 


The  elastic  properties  of  the  armature/rail  system  need  to  be  considered.  The  lateral 
force  on  the  rails  by  the  magnetic  field  is  a  load  that  moves  with  the  velocity  V  of  the 
armature.  It  is  known^  that  moving  loads  on  elastically  supported  beams  can  excite 
disturbances  in  the  beam  and  more  important,  that  there  is  a  critical  velocity  V *  at  which 
the  energy  does  not  radiate  away  from  the  moving  load,  but  builds  up  in  its  vicinity.  In 
the  absence  of  damping,  the  response  becomes  infinite  as  V  approaches  V*.  The  value  of 
this  critical  velocity  is 


l/* 


(AkEI\ t 

WA>) 


(27) 


where  E  is  Young’s  modulus,  /  is  the  moment  of  inertia,  p  is  the  density,  A  is  the  cross- 
section  area  and  k  is  the  restoring  force  per  unit  displacement  of  the  beam  in  the  lateral 
direction.  When  V  reaches  V* ,  the  displacement  of  the  rail  may  be  sufficient  to  break  the 
contact  between  the  armature  and  rail  and  cause  transition.  Research  on  this  possibility 
will  be  reported  later. 
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Comments  on  the  Paper 


This  report  embodies  the  paper  that  I  am  submitting  to  the  8th  EML  Symposium 
to  be  held  in  April  in  Baltimore.  It  is  attached  below  and  may  undergo  minor  revisions 
before  presentation.  It  is  based  on  the  theory  presented  in  my  first  Progress  Report,  which 
gave  an  account  of  a  melt-wave  model  for  transition  in  solid-armature,  rail-gun  launchers. 
This  report  contained  a  number  of  assumptions  concerning  the  electrical,  mechanical  and 
thermal  interactions  between  the  rail  and  the  armature.  These  were  (i)  perfect  electrical 
contact,  (ii)  the  removal  of  liquid  metal  by  viscosity,  (iii)  in  addition  to  the  (obvious) 
boundary  condition,  j*  =  j*  ,  on  the  normal  component  of  the  current  density,  a  boundary 
condition  on  the  transverse  current  density  at  the  armature/rail  interface  of  the  form 
VrJx  =  Vajxt  where  the  subscripts  ‘a’  and  V  denote  armature  and  rail  and  77  is  the 
resistivity,  (iv)  adiabatic  conditions  except  possibly  in  the  region  upstream  of  the  melt- 
wave,  (v)  that  the  width  of  the  gap,  h  between  the  armature  and  the  rail  is  based  on  an 
electrical  scale-length  and  (vi)  (dependent  on  (iii))  that  this  scale  length  is  8a  =  ria/(p,0V). 

In  preparing  the  EML  paper,  I  decided  to  replace  the  dimensional  arguments  used  to 
obtain  (vi)  by  a  mathematical  solution  of  the  equations  satisfied  by  the  current  density. 
The  dimensional  argument  was  based  on  the  two  boundary  conditions  in  (iii),  which  lead 
naturally  to  two  distinct  length  scales  for  change  in  the  currents  within  the  armature, 
namely  8r  and  6a.  The  thickness  h  was  then  identified  with  the  length  scale  of  the  normal 
current,  jy  in  the  OT-direction,  which  proved  to  be  6a.  But  this  is  not  support  by  the 
mathematical  solution,  which  shows  that  there  can  be  only  one  length  scale  within  the 
armature  (see  equation  (9)  below).  This  length  scale  is  that  imposed  on  the  armature 
currents  by  an  outer  boundary  condition  in  the  rail;  it  is  therefore  6r,  not  6a.  The  basic 
error  in  the  first  report  was  to  follow  Parks  [3]  and  adopt  the  second  boundary  condition 
in  (iii). 

The  second  boundary  condition  followed  from  the  transverse  component  of  Ohm’s 
law,  r) j  =  E  +  vxB,  applied  to  each  side  of  the  interface,  and  the  steady  state  equation 
VxE  =  0,  integrated  across  the  boundary,  which  reputedly  shows  that  the  transverse 
component  of  the  electrical  field  is  continuous.  At  first  sight  this  is  plausible,  but  it  is 
not  valid  for  a  reason  I  gave  in  my  book,  Principles  of  Magnetoplasma  Dynamics ,  OUP, 
1987,  (pp.  116-17)  and  had  overlooked.  Boundary  conditions  are  obtained  by  applying 
transition  rules  to  the  field  equations  and  then  applying  physical  conditions  for  setting 
surface  intensities  equal  to  zero.  By  surface  intensities  is  meant  delta  function  infinities 
in  the  property  in  question,  a  familiar  example  being  the  current  sheet  j*  appearing  in 
the  transition  rule  nx[B]  =  ^0j*,  obtained  by  integrating  VxB  =  /i0j  across  a  surface 
with  unit  normal  n.  Only  if  we  can  adduce  physical  reasons  for  j*  being  zero,  such  as 
a  significant  resistivity,  can  we  infer  the  boundary  condition  nx[B]  =  0.  Otherwise  the 
condition  merely  gives  us  a  means  of  determining  the  strength  of  the  current  sheet — it  is 
not  a  boundary  condition. 

When  the  field  equation  VxE  =  0  is  integrated  across  a  surface  layer,  the  result  is 

nX[E]  +  V|,  XE*  =  0, 

where  V||  is  an  gradient  operator  lying  in  the  surface  and  E*  is  an  electric  field  surface 
intensity.  The  physical  origin  of  E*  is  an  electric  dipole  layer  of  strength  r*  say,  which 
gives  E  =  nr* .  If  the  length  scale  of  interest  in  the  physical  problem,  £  say,  is  greater 
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than  a  Debye  length,  it  is  not  possible  to  set  r*  equal  to  zero.  The  Debye  length  is  very 
small  on  an  MHD  scale  of  length  so  in  rail-gun  physics  we  cannot  set  E*  equal  to  zero, 
which  in  turn  means  that  we  are  unable  to  assume  that  nx[E]  =  0.  It  is  therefore  not 
true  that  rjaj%  =  r\rjTx  and  therefore  the  paper  has  only  one  length  scale  for  the  armature 
current,  namely  6r. 

In  my  second  report  ( The  Contact-Spot  Model  of  Transition)  it  was  explained  that 
the  liquid  model — i.e.  the  model  in  which  most  of  the  melted  metal  remained  filling  the 
gap — is  valid  only  if  the  rate  at  which  the  volume  of  the  pool  is  increased  by  melting 
exceeds  the  rate  at  which  it  is  removed  in  a  boundary  layer  by  friction  with  the  rail.  The 
model  was  rejected  since  the  condition  on  the  rates  is  not  met.  With  the  entrainment 
model,  the  fluid  is  never  ‘deep’,  so  the  ‘boundary  layer’  is  not  a  relevant  concept.  It  seems 
very  reasonable  to  expect  that,  as  fast  as  the  melt  is  created  by  ohmic  heating,  it  will 
be  removed  simply  by  adherence  on  the  rail.  As  shown  at  the  end  of  the  attached  EML 
paper,  this  model  gives  a  thickness  Sm  of  solidified  Aluminum  on  the  rail  of  about  20  ^m, 
a  value  that  is  supported  by  measurements.  I  have  included  one  other  important  change 
at  this  point— I  have  allowed  for  pressure  between  the  armature  and  rail  by  admitting 
repeated  traverses  of  the  armature  by  the  melt  wave.  Transition  is  defined  as  occurring 
when  removal  of  the  armature  surface  reduces  the  pressure  to  zero. 

Returning  to  the  list  in  the  first  paragraph  above,  I  have  removed  the  adiabatic  as¬ 
sumption  of  (iv),  a  change  that  generates  the  ‘G’  term  in  equations  (18)  and  (20)  of  the 
paper  below.  It  proves  to  be  only  a  small  correction  in  (18),  but  more  significant  in  the 
expression  (20)  for  the  velocity  of  the  armature  at  transition. 

A  more  important  matter  is  the  question  of  what  process  determines  the  thickness, 
h,  of  the  gap  between  the  solid  armature  and  the  rail  behind  the  advancing  melt- wave. 
In  my  first  report  I  adopted  (v),  i.e.  accepted  an  electrical  scale  length  for  h,  but  in  my 
second  report  (mainly  concerned  with  the  nature  of  the  electrical  contact)  I  changed  to  a 
thermal  scale  length  by  an  uncertain  argument  that  no  longer  convinces  me.  It  may  be 
correct,  but  in  the  EML  paper,  I  have  thought  it  prudent  to  return  to  the  electrical  scale 
length.  Supporting  this,  there  is  the  fact  that  the  time-scale  for  electrical  changes  over 
a  given  length,  being  inversely  proportional  to  the  magnetic  diffusivity,  is  much  smaller 
that  the  time  scale  for  thermal  changes  over  the  same  distance.  I  expect  to  return  to  this 
question  in  later  reports. 
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The  Current  Melt- Wave  Model 


L.  C.  Woods 

Institute  for  Advanced  Technology,  University  of  Texas  at  Austin 
&  Mathematical  Institute,  University  of  Oxford 


Abstract  It  is  generally  accepted  that  the  velocity  skin  effect  concentrates  the  cur¬ 
rent  crossing  the  armature/rail  interface  at  the  rear  corners  of  a  solid  armature.  Ohmic 
heating  melts  these  corners  and  the  molten  metal  is  removed  by  being  entrained  on  the 
rail.  As  the  resistivity  of  the  resulting  gap  exceeds  that  of  the  solid  phase,  the  current  is 
deflected  forwards  into  the  solid  region.  The  mechanism  results  in  a  ‘melt-wave’  moving 
from  the  back  to  the  front  of  the  armature,  leaving  a  high  resistance  plasma  gap  of  width 
h,  between  the  armature  and  the  rails.  Transition  to  a  high  contact  voltage  is  therefore 
a  likely  outcome  when  the  melt  wave  reaches  the  front  of  the  armature,  although  if  the 
armature  is  compressed  by  the  rails,  several  successive  meltings  along  the  armature  sur¬ 
face  will  be  required  to  reach  transition.  In  this  paper  we  give  a  theory  for  the  speed  vm 
at  which  the  melt  wave  advances  relative  to  the  armature.  Our  model  differs  from  those 
proposed  by  other  authors  in  that  h  depends  only  on  the  electrical  skin  depth  in  the  rail 
and  since  vm  oc  h-1,  our  value  for  vm  is  different.  An  expression  for  the  armature  speed 
Vc  at  which  transition  can  occur  is  also  derived. 


1.  INTRODUCTION 

The  main  assumptions  on  which  our  model  depends  are  (i)  perfect  electrical  contact 
between  the  armature  and  rail  ahead  of  the  melt  wave  and  (ii)  the  rapid  removal  of  the 
liquid  metal  by  entrainment  on  the  rail.  The  first  assumption  is  not  very  satisfactory, 
since  the  electric  current  is  actually  transported  across  the  gap  between  rail  and  armature 
through  small  contact  spots  [1]  [2],  but  the  model  will  at  least  serve  to  establish  an  ideal 
value  for  the  time  to  transition.  Also,  if  the  pressure  between  the  surfaces  is  ‘large’,  the 
spots  will  increase  in  size  and  number  under  this  normal  force  and  the  electrical  contact 
will  begin  to  resemble  the  ideal  case. 


Figure  1:  Melt- wave  model 
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The  melt-wave  model  is  illustrated  in  Fig.  1.  The  wave  is  advancing  with  a  speed 
vm  relative  to  the  armature,  moving  from  D  to  B,  and  when  it  reaches  f?,  provided  the 
gap  h  remains  large  enough  not  to  be  suppressed  by  the  normal  pressure  between  the 
surfaces,  transition  will  result. 


An  expression  for  vm  can  be  deduced  from  the  energy  equation  applied  to  the  solid 
region  of  the  armature  extending  in  front  of  the  melt  wave.  We  shall  show  that  it  is 
given  by 


Vm  — 


afuioPV  2apoGV 

2r)TpQwi  irrir 
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2ot(j,oV  ^ 
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1.148)  , 


(1) 


where  Q  =  cp{Tm  -T0)+Lm,  G  =  ^lXa  [l  +  (~)5]  >  (2) 


and  for  completeness,  we  have  included  the  expression  for  h.  In  these  expressions  rja 
is  an  average  resistivity  in  the  solid  armature,  7]r  is  the  rail  resistivity,  I  is  the  total 
current  flowing  through  the  armature,  p o  is  the  permeability  of  free  space,  p  is  the 
armature  density,  w  is  the  depth  of  the  surface  of  the  armature  in  contact  with  a  rail, 
Tm  is  the  melting  temperature,  To  is  the  laboratory  temperature,  Lm  is  the  enthalpy  of 
fusion  (latent  heat),  cv  is  the  specific  heat  at  constant  pressure,  Xa?  Xr  are  the  thermal 
diffusivities  in  the  armature  and  rail  and  V  is  the  speed  of  the  armature.  The  armature 
speed  Vc  at  which  transition  occurs  is  given  in  equation  (20). 

While  Parks  [3]  gives  an  accurate  description  of  the  phenomenon,  he  obtains  values 
for  the  scale  length  for  changes  in  the  magnetic  field  that  appear  to  be  in  error.  His 
expression  for  the  melt-wave  gap,  h,  includes  an  additional  factor  0.45(rya/r/r)2,  and 
since  vm  oc  h-1,  his  formula  for  vm  disagrees  with  ours.  Also  the  conflict  between  his 
figures  1  and  2  is  puzzling,  as  is  his  remark  that  the  presence  of  the  gap  between  the 
armature  and  the  rails  will  not  appreciably  alter  the  distribution  of  the  current. 

James  [4]  adopts  a  model  in  which  the  melted  metal  from  the  armature  is  only  slowly 
removed  by  the  rail.  It  therefore  fills  the  gap  between  the  armature  and  rail,  which 
allows  some  current  to  flow  across  the  gap.  To  some  extent  this  weakens  the  melting 
process  and  it  also  means  that  melting  alone  is  not  sufficient  to  yield  transition.  James 
therefore  defines  the  current  wave  to  be  a  vaporisation  wave.  To  modify  equation  (1)  to 
apply  to  this  case,  we  need  to  increase  the  value  of  Q  to  include  the  additional  internal 
energy  between  the  melting  and  vaporisation  temperatures  plus  the  heat  of  vaporisation. 

Barber  and  Challita  [2]  have  no  confidence  in  melt- wave  morels,  which  they  find 
irrelevant  because  the  assumption  of  a  perfect  contact  surface  between  rail  and  armature 
is  never  attained  in  practice.  The  electric  current  flows  through  discrete  spots  that 
occupy  only  a  small  fraction  of  the  surface,  and,  they  argue,  this  considerably  weakens 
the  influence  of  the  velocity  skin  effect.  Their  treatment  of  armature  contact  transition, 
is  based  on  ‘adiabatic  accumulation  of  resistive  dissipation  in  contact  spots’.  Whether 
the  contact  spots  are  really  adiabatic  depends  very  much  on  their  size,  which  is  poorly 
known.  We  agree  with  some  features  of  their  model.  Ours  is  an  idealized  description  with 
which  to  compare  experiments  and  more  physical  representations  of  electrical  contact. 
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2.  THE  CURRENT  DENSITY 


We  shall  adopt  a  two-dimensional  model,  as  depicted  in  Fig.  1,  with  axes  fixed  in 
the  armature.  The  magnetic  field  B  is  directed  along  OZ,  the  rails  move  parallel  to  OX 
with  speed  V  and  the  electric  current  j  lies  in  the  (X,  Y)-plane.  Unit  vectors  along  the 
axes  are  denoted  by  x,  y,  z.  The  equation  determining  the  evolution  of  the  magnetic 
field  in  an  incompressible  medium  is  [5] 

dB 

-gT  +  Ve-  VB  =B-  Vve-KV2B,  (3) 


where  ^  =  rj/y,0  is  the  magnetic  diffusivity  and  ve  is  the  electron  fluid  velocity, 

V  x  —  j/ene  (rail) 

—j/ene  (armature) , 


in  which  ene  is  the  electron  charge 

density. 

With  the  relations 

dBz 
dx  ~ 

dBz 

dy  -Mz, 

(5) 

we  obtain 

ve  •  VB  = 

J  -yojyVz 

(rail) 

\  Oz 

(armature) . 

We  shall  assume  that  in  a  frame  fixed  in  the  melt-wave  front,  the  thermal  and  electric 
fields  are  in  a  quasi-steady  state,  by  which  we  mean  that  on  a  time  scale  short  compared 
with  the  transit  time  of  either  the  armature  from  breech  to  muzzle  or  the  melt  wave 
through  the  armature,  the  fields  will  appear  static.  This  approximation  will  evidently 
fail  near  either  end  of  these  transits.  In  a  typical  experiment  at  IAT,  Austin,  we  find 
that  the  diffusion  time  for  the  electric  field  through  a  j  inch  copper  rail  is  about  1  ms, 
which  is  not  much  less  than  the  armature  transit  time.  So  a  genuine  steady  state  in  the 
electric  field  is  never  attained. 

Since  B  •  V  =  0,  in  steady  conditions  (3)  yields 

,  (dL  +  jL\R  _  J  -*>3vv  MO 

V^2  dy2)  2  |  0  (armature)  _ 


We  shall  use  superscripts  V  and  ‘a’  to  denote  rail  and  armature  values.  By  (5)  we 
can  write  the  equations  as 


dj, 


+  v^  =  0’ 


We  also  have  V  •  j  =  0  in  both  domains,  i.e. 


dj%  dJy 

dy  dx 


=  0. 


dx  dy  ’ 


dj%  djy 
—  +  —2-  =  0 
dx  dy 


(6) 

(7) 
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Eliminating^  from  (6)1  and  (7) j ,  we  get 


where  6r  =  £r/V.  Notice  that  with  x'  =  x/6r,  y'  =  y/6r,  the  equation  takes  a  dimen¬ 
sionless  form. 

We  now  assume  that  the  wave  front  lies  on  x  =  0  and  that  the  rail/armature  interface 
lies  on  y  =  0.  As  we  are  interested  in  a  solution  near  the  origin,  the  geometry  suggests 
that  we  solve  (8)  by  a  separation  of  variables,  i.e.  by  writing  jr  =  X(x)Y (y).  We  require 
jy  to  vanish  at  x  —  oo.  The  solution  is  therefore  the  sum  over  the  separation  parameter 
m  of  terms  of  the  the  type 

fy  =  Am  cos (my/Sr)eax/6r ,  a  =  1{(1  +  4m2)s  -  1}  (x  <0,  y  <  0).  (9) 

To  determine  m  we  consider  the  diffusion  of  the  electric  field  into  the  rail  as  viewed 
from  the  armature.  In  a  time  interval  6t,  a  current  element  diffuses  a  short  distance 
y  =  -(irfit)?  into  the  rail  as  it  is  convected  a  distance  x  —  VSt  along  OX.  The  element 
therefore  moves  at  a  small  angle  6  =  (£rSt)? / (V6t)  to  OX.  Hence  its  displacement  along 
OY,  as  seen  from  the  armature,  is  0(£r6t)?  =  £r/V  —  <5r.  It  follows  that  jy  must  fall  to 
zero  at  y  =  -Sr.  Hence  by  (9),  m  =  7r/2  and  a  =  1{(1  +  n2)?  -  1}  «  1.148. 

On  integrating  V  •  j  =  0  across  the  interface  y  =  0,  we  obtain  the  boundary  condition 
Jy(x.  0)  =  jy(x,  0),  whence  from  (9) 

jy(x,  0)  =  j0eax^  (-oo  <x  <  0),  (10) 

where  j0  is  the  current  density  at  the  origin.  The  appropriate  continuation  of  this 
boundary  value  into  the  armature  is 


y ) 

jy(x,  y ) 


sin  (ay /6r) 
cos(ay/ST) 


joe 


ax/6r 


(— oo  <  x  <  0,  0  <  y  <  h), 


(11) 


which  satisfies  (6)2  and  (7)2.  We  chose  h  to  be  the  distance  at  which  the  current  flowing 
into  the  armature  first  falls  to  zero,  i.e.  where  ah/6r  =  ir/2.  Hence 

7C& 

h  =  «  1.368  6r  =  \kI,  (l  =  Sr/a) .  (12) 


From  (11) 

jy  (z,  0)  =  jo e®7*,  j2(x,  y)  =  j%(?x/e ,  (13) 


where  j  is  the  armature  current  density.  The  total  current  passing  through  the  surface 
y  =  0  is  therefore 


fo 

I  =  W  /  jydx  =  wtjo,  (14) 

J  —  OO 


where  w  is  the  depth  of  the  surface  of  the  armature  in  contact  with  the  rail. 
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3.  THE  SPEED  OF  THE  MELT  WAVE 


We  start  with  the  energy  equation 

pepDT  kV  T  +  rjj  pLrnvm6(xm),  (15) 

where  D  =  d/dt  +  v  •  V  is  the  convective  time  derivative,  k  is  the  thermal  conductivity, 
6{xm)  is  the  delta  function  and  the  other  variables  have  been  defined  in  the  Introduction. 
Our  assumption  that,  in  a  frame  fixed  in  the  wave  front,  conditions  are  steady,  allows  us 
to  write  DT  =  vmdT/dx.  Equation  (15)  now  takes  the  form 

dT  ( d2T  d2T\ 

pCpVm~dx  =  Ka  (  'dx2  +  ~dy2)  +  Vajo^*  ~  pLmvm6(0) .  (16) 


This  equation  is  next  integrated  over  -oo  <  x  <  0  and  0  <  y  <  h.  As  dT/dx  is  zero 
at  x  =  0  and  x  =  —  oo,  the  thermal  conductivity  term  yields 


f  «. 

J—  OO 


i  h 


dT 

L5yJ0 


dx 


rO 

-KaiK  +  K  )  /  (T-  To)  dx  «  -«0(1  +  Aa/Ar)(Tm  -  T0), 

J  —oo 


where  Aa,  Ar  are  the  thermal  skin  depths,  (xato)$  and  (xf<o)*-  With  the  help  of  (12) 
through  (14),  the  integrated  equation  yields  the  solution  given  in  (1)  and  (2),  viz. 


Q7?q^qJ2F  _  2ap0GV 
2t]t  pQw2  npr 


(.  T^h 

V  2  apoV' 


a 


(17) 


where  Q  =  cp(Tm  -  T0)  +  Lm  ,  G  =  7o)ya[l  +  (^^]  .  (i8) 


The  acceleration  of  the  armature/projectile  combination  is  given  by  the  well-known 
formula 


dV 

dt 


L'l 2 
2  M 


-T, 


(19) 


where  M  is  the  total  mass  being  accelerated,  V  is  the  induction  gradient  and  T  is  the 
sum  of  the  various  drag  forces  acting  on  the  combination.  We  shall  assume  that  the 
armature  is  under  pressure  from  the  rails  and  that  until  a  width  /  h  is  removed  from 
each  side  of  the  armature,  electrical  contact  is  maintained.  Eliminating  I  between  (17) 
and  (19)  and  integrating  vm  =  dxm/dt  over  the  time  tm  =  fa/vm  that  it  takes  the  melt 
wave  to  travel  the  length  a  of  the  armature  f  times  and  rearranging,  we  get 


V2 

v  c 


2r}rpQw2Li 

arjafj,0M 


(fa  +  -  2  [tm  VT dt , 

V  7T?7r  /  Jo 


(20) 


where  Xm  is  the  distance  between  the  breech  and  the  armature  when  the  melt  wave  has 
destroyed  the  electrical  contact.  We  identify  Vc  as  the  armature  speed  when  transition 
occurs. 
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The  thickness  6m  of  the  melted  armature  metal  deposited  in  the  rail,  can  be  calculated 
as  follows.  On  the  assumption  that  the  mass  entrained  on  the  rail  equals  the  mass  melted 
in  the  wave  front,  we  obtain 

Ps'VrrM  ~  Pi(V  +  Vm)fim  i 

where  the  subscripts  s  and  l  denote  the  solid  and  liquid  phases.  Since  vm  V,  this 
yields 

sm  =  ^  h.  (21) 

v  Pi 


4.  APPLICATION 

Suppose  that  the  armature  metal  is  Aluminium  and  the  rails  are  Copper  then  we 
need  the  following  physical  values,  all  expressed  in  SI  units: 

rja  =  7.2xl0~8,  i ir  =  1.7xi(r8,  ps  =  2660,  Pl  =  2410,  Tm  -  T0  =  635, 

Q  =  9.67xl05,  cp{Tm-T0)/Q  =  0.589,  Xa  =  1.07xl0"4,  Xr  =  1.16xl0-4. 

With  these  values,  equations  (17)  and  (20)  become 

vm  =  1.19xl0~15^  -6.68x10-3P,  h  =  1-85x10  2 
wz  V 

and  Vc2  =  8.42xl014^^(/a  +  6.67xl0-3Xm).  (23) 

where  we  have  omitted  the  drag  term.  The  thickness  of  the  deposited  Aluminium  follows 
from  (21): 

6m  =  1-l^h.  (24) 

To  illustrate  the  theory,  we  shall  apply  it  to  a  rail-gun  experiment  carried  out  in 
the  Institute  for  Advanced  Technology  [6].  The  armature’s  initial  (free)  width  of  26  mm 
is  compressed  to  24.8  mm  to  fit  between  the  rails.  It  follows  that  /  =  1.2xl0_3/2h. 

Other  typical  values  are:  w  =  12.7  mm,  a  =  20.3  mm,  I  —  2.5x10s  A,  V  = 

500ms-1,  and  L'/M  =  2.2xl0~5.  These  yield: 

vm  =  0.45  V  =  225 ms-1,  h  =  37 pm,  Vc  =  l.OOxlO3  ms-1, 

6m  —  18 /im,  tm  =  1.46  ms,  Xm  =  0.56  m, 
which  values  are  in  broad  agreement  with  experiment. 
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An  Armature  /  Rail  Instability 
Summary 

It  is  known  that  the  movement  of  a  load  along  an  beam,  supported  on  an 
elastic  foundation  will  excite  the  beam  into  transverse  vibrations  that  are  par 
ar  y  enhanced  when  the  speed  V  of  the  load  is  close  to  the  velocitv  of  free 
waves  in  the  beam  The  lowest  velocity  at  which  a  free  wave  can  propaga  e  defies 
the  critical’  velocity,  which  is  given  by  p  8  etmes 

C  l  P2  )  ’  (1) 

where  E  is  the  modulus  of  elasticity  of  the  beam  r  ;Q  * K  > 

Z t  f  “d  f '  f  is;he  ^  >-^t“rc°ef 

is  damping  ofT^l 

disturbance  increases  without  limit,  u  we  have  an'  ^“n^cf  ““ 

Experiments  with  rail  guns,  in  which  an  Aluminium  armature  moves  along 
copper  rails  at  a  speed  exceeding  a  critical  value  V*,  have  revealed  a  tvnn  nf  u 
damage  known  as  ‘gouging'.  Following  the  fast  eoug™e t a  luenle 
gouges  on  each  rail  that  appear  to  be  periodic,  with  a  wave  length  that  h  t 
ens  with  increasing  distance  along  the  rail.  The  magnetic  force  propelling  the 
armature  has  components  in  the  direction  of  the  rails.  Thus  the  armature  k 

1  I  Thm°T8  °Ver  *W0  “PPer  b6ams'  each  s“PP°rt<=d  on  a  foundatfon  mate 
his  description  raises  the  possibility  of  a  connection  between  V  for  the 

ra,l/foundat,on  system  and  the  speed  V  at  which  gouging  commences.  ‘ 

The  amplitude  of  waves  propagating  along  the  rails  at  speeds  of  V  —  OV 
depends  on  the  physical  parameters  defined  above  and  on  thT  ratfa  fl  The  " 
are  features  of  these  waves  in  the  neighbourhood  of  9  ~  1  that  ■  Tbere 

significant.  We  shall  describe  them  for  the  ci  of  7  Partlcular‘y 

conclusions  are  unaltered  for  a  small  amount  nf  A  •  amped  system— the 

the  load,  the  waves  appear  itfa  g  “Ttfal"8  T 

position  of  the  load.  Let  *  =  0  define  the  position  of  the  load,  wh,ch  is  burned 

acting  on  the  beam.  In  the  neighbourhood  of  6  =  1,  '°ad  ? 


9  <  1  : 


cos  krx 


2V2F(\-Q2))i  ' 


9>  1  : 


sin  krx 


2 F  ( 9 4  -  l)* 


where  kc  =  ( F/EI)i . 

In  interpreting  (2),  we  need  to  remember  that  there  are  two  rail,  •  * 

beams,  and  that  the  point  loads  correspond  to  the  resultant  f  ’  ^  tW° 

surfaces  of  the  armature  in  contact  with  the  rail  For  In  AluiW  ^  ^ 

we  shall  assume  that  the  rails  are  displaced  symmetrically  abou”  "he”  Zhnel’ 
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of  t!1.e  ®un’  In  thls  case  the  armature  moves  steadilv  along  L.  without  transverse 
oscillations.  (The  case  of  antisymmetrical  displacements  will  be  discussed  at  the 

the  h*  115  rCP°(  ,For  subcntIcal  speeds’  W  <  1).  the  first  of  (2)  shows  that 
he  displacement  is  in  the  same  direction  as  the  load  on  each  rail  Hence  with 

symmetncal  displacements  (see  Fig.  2),  the  load  tends  to  move  the  rails  ap« 
and  allow  freer  passage  for  the  armature  along  the  gun  Provided  this  J  P 
is  not  so  great  that  it  interrupts  the  flow  of  current' acroj,  ,L  ^  <?anS10n 
causes  transition,  there  will  be"  no  adverse  lsequeli  ' „  tt  7 

separate  the  load  is  reduced  until  a  balance  is  attained  between  the  force  o„”he 
rails  and  the  response  of  the  elastic  foundations  iP  rce  on  the 

stable.  However  for  speeds  just  supercritical,  the  second  o ^  “als"  th^f " 
x  >  0’  a  position  just  in  front  of  the  resultant  force  P  th0  •,  eais  that  for 
inwards.  Hence  P  is  both  increased  and  moved  forwards’  4s  the  fo^  SplaCed 
the  inwards  displacement  becomes  greater  and  P  ic  inn  At  .  f°rCe  lncreases» 
the  force  between  the  rails  and  arm  1“  “Ifa «  The  T  T  ^  ^ 

softened  by  the  ZnTZtZZ  ZZreZ  ThT  t  1“*"  ba 

phenomenon  initiates  the  periodic  gouging  of  the  rails  ie  thaM/' that  th‘S 

to  be  plausible.  The  actual  details  of  the  /  •  ’  that  Vc  =  V  ’  aPPears 

in  this  report  we  are  con 

«  ^  VZKe  77 is 

a~Sder»dl^^^ 

amplitude  may  result  in  shallower  gouge  marks  but  both  reduCtl0f  ln  wave 

depth  of  gouging  require  a  more  detailed  mechanism  for  ^  T  and  the 
we  are  presenting  here.  Sm  f°r  their  exPlanation  than 

One  of  the  main  difficulties  in  applying  the  thpnrv  Hoc  •  a 
response  of  the  foundation,  represented  by  F  in  m  The  ?  d<®Tmmin«  the 
to  treat  the  foundation  as  being  a  seoumcc  ,,f  ,  '  ,~h  al  aPPr»nch  is 

is  evidently  a  poor  model  of  the  aS  T  1,ndePendc”t  'springs',  but  this 
continuous  foundation  material  The  mu'  P  . pr0Cess  of  the  resPonse  of  a 
a  method  of  determilg  “  To  d^hisT TTl  C°"tribution  in  this  report  is 
of  a  displacement  that  is  distributed  over  thTwhltf “ 
approximations  are  unavoidable  For  the  tat  •  th6  toundatlon  material; 
yields  V*  =  1380 m/s,  in  ^ 

for  ^‘undX  ^  f°Und  ^  value 

Copper  armature  is  3.32  times  heavier  than  an^T  — °S  &t  ab°Ut  600  m/s-  The 

shape  and  this  appears  to  be  the  only  difference  fo  thrth™^11"6,^  ^  ^ 
the  mass  of  the  armature  has  no  role  t  ™  theory  outlined  above, 

of  the  rail  are  induced  by  the  greater  mas-?  *  °'nttsyTnmetric  displacements 
We  shall  discuss  this  in  §5  §  maSS’  &  mUCh  loWer  criticaI  aPeed  results. 
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1.  The  Critical  Armature  Velocity 


We  base  our  account  on  the  paper  by  J.T.  Kenney  [1],  who  cites  several  earlier 
works.  A  useful  survey  of  the  theory  is  presented  by  Graff  [2],  There  follows  an 
outline  of  those  parts  of  Kenney:s  paper,  relevant  to  rail  guns. 

Let  x  denote  the  distance  along  the  beam,  y  denote  its  transverse  displacement 
and  w  be  the  load  per  unit  length,  as  illustrated  in  Fig.  1.  If  M  is  the  bendimr 
moment,  then  the  classical  relations  ° 


il  (EIdV 

dx 2  \  dx2 


a2M 

'  dx2 


=  w, 


(3) 


apply,  provided  the  deflection  y  is  relatively  small.  The  several  loads  are  given 

w  =  P(x,  t)  —  p^-^-  —  c—  —  Fv 
V  Pdt2  dt  V ’ 

where  P(x ,  t)  is  the  applied  force  and  the  other  right  hand  terms  are  the  inertia 
ampmg  and  foundation  loads.  We  have  introduced  p,  the  mass  per  unit  length’ 
e  damping  coefficient  c  and  the  spring  constant  per  unit  length.  F.  The  as- 

bvTefor^V8^  the  reSP°nSe  °f  the  founda-tion  can  be  represented  locally 
by  the  force  F.  This  is  serious  approximation,  to  which  we  shall  return  in  fj3 
The  differential  equation  for  the  beam  is  therefore 


-d2y  dy 

EIdx 4  Pdt2+Cdi  +  Fy  =  p(x ’ 


(4) 


figure  1:  Deflection  of  beam  by  a  load 

alnnWeiKarVntereStld  ^  a,  possibIe  resonance  between  an  applied  force  moving 
along  the  beam  with  a  velocity  V,  i.e.  a  force  of  the  type  P(x  t)  =  P(r  -  Vt\ 

and  ‘natural’  waves  that  could  exist  in  the  beam  in  the  absence  of  an  applied 
orca  Of  course  the  force  is  required  to  excite  these  waves,  but  if  in  add-on 
this  force  is  resonant  with  natural  frequencies,  we  can  expect  large  amplitudes 
to  result.  To  determine  the  natural  waves,  we  seek  wave-type  solutions  of  the 
homogeneous,  undamped  form  of  (4),  viz.  ™ 


with 


y  =  Y exp {i(kx  -  tut)}  (Y  =  const.). 


The  dispersion  relation  is 


Elk4  -  pco2  +  F  =  0 . 

Define  a  =  ( EI/F)k 4,  0  =  (p2 /FEI)(u/k)4  , 


(6) 

(7) 


then  the  wave  length  is  A  =  2? r/k  =  2tt (El/Fa)* 

and  (6)  becomes  a  —  (a/?)i  +  l  =  o , 


i.e. 


P  —  OL  +  2  +  a 


(9) 


u 


Since  the  phase  velocity  of  the  waves  is  vp  =  u/k,  we  see  that  0  determines  ^ 
To  ensure  propagating  rather  than  evanescent  waves,  we  require  k  to  be  real  and 
therefore  a  to  be  positive.  The  lowest  phase  velocity  occurs  when  the  right-hand 
side  of  (9)  is  a  minimum,  which  occurs  at  a  =  1,  where  0  =  4.  It  follows  that  the 
lowest  velocity  at  which  a  free  wave  can  propagate  is  vp  =  Vc,  where 


V, 


-  (4FEI\* 

~\P2) 


(10) 


Therefore  resonance  occurs  when  the  load  moves  with  this  speed  and  from  (7) 
and  (8),  the  corresponding  frequency  and  wave  length  are 

u>c  —  (2T1/ p)1 ,  Ac  =  2'K/kc  =  2ir(EI / F)^ ,  kc  =  (F/EI )i  .  (11) 


2.  Free  standing  waves 


The  next  step  in  the  analysis  is  to  transform  into  a  coordinate  system  moving  with 
the  load  and  then  to  investigate  standing  waves  in  this  frame.  This  is  effected 
by  writing  *!=**-  Vt,  Vl  =  y(x  -  Vt)  =  yi(Xl)  and  then  setting  all  partial 
derivatives  with  respect  to  time  equal  to  zero.  This  simplifies  equation  (5)  to 

Ely™  +  pV2y"  +  Fy  =  0  ,  (i2) 

where  the  primes  denote  differentiation  with  respect  to  xx.  The  load  P  does  not 
appear  on  the  right  hand  side,  since  it  is  now  concentrated  at  the  origin  and 
enters  the  solution  as  a  boundary  condition,  there  being  a  discontinuity  in  shear 
across  the  origin  of  magnitude  -P/EI.  (Of  course  we  could  have  included  it  in 
(12)  as  a  delta  function  at  the  origin.)  When  damping  is  included,  the  general 
solution  is  complicated  algebraically  (see  equation  (14)  of  [1]),  but  for  our  present 
purposes  it  is  sufficient  to  consider  the  undamped  case. 
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We  adopt  the  definitions 

e  =  V/VCy  G±  =  ^{(^  +  I)i±(^_i)i}>  k  =  ^(02  +  1)^  (13) 

then  the  solution  of  (12)  can  be  expressed  (cf.  equations  (14a)  and  (14b)  of  [1]) 

9<l-  V  =  ^^;^p{±ks(l-e2)^}  + 

2V2F  l  y/2. K  (l-02U/> 


_  _Pk<  sin jkcG^x) 
J  ”  2 F  (04-l)*G± 


9>  1:  =  _£^c  sm[KcU~x)  x<0.  _sign 

2F  (0<-\)hG±  x  >  0;  +sign  ’  (14) 

where  the  displacement  y  of  the  beam  is  positive  if  in  the  same  direction  as  the 
force  applied  to  the  beam  and  we  have  introduced  ks  =  kc/V 2,  the  wave  number 
for  the  static  case,  i.e.  the  wave  number  at  6  =  0. 

°f  thiS  S0luti0"  is  not  valid  6  =  1.  Damping  in  the  response  of 

the  foundation  is  required  to  suppress  the  infinity.  Also  observe  that  y  and  all 
its  derivatives  are  continuous  at  the  origin.  In  the  neighbourhood  of  0  =  1  „e 
obtain  the  results  given  in  (2).  ’  e 


rail  distorsion  at  the  critical  speed 


Figure  2:  Armature  gouging  rails 

The  difference  between  the  type  of  displacement  for  6  <  1  and  9  >  1  is  the  kev 
feature  so  far  as  gouging  is  concerned.  For  subcritical  speeds,  the  displacement 

smVSierr  ^  '°ad  "  SUffiCiently  redU“d  for  “  equilibrium 
state  to  be  achieved.  The  rail/armature  system  is  stable.  However  for  speeds 

just  supercritical  for  *  >  0,  the  rails  are  displaced  inwards.  Hence  the  force  i 

both  increased  and  moved  forwards.  As  the  resultant  force  increases,  the  inwards 

displacement  becomes  greater  and  the  force  is  further  increased.  The  electrical 

contac  between  the  rail  and  armature  will  have  least  resistance  where  the  force 

greatest  so  we  can  expect  ohmic  heating  to  have  softened  the  copper  at  this 

point.  (At  a  temperature  of  400°  C,  copper  has  only  about  a  fifth  of  its  room 
temperature  tensile  strength.)  00m 

Gouging  may  occur,  as  indicated  in  Fig.  2.  Viewed  in  a  frame  with  thp 
armature  stationary,  the  rail  will  'snake'  by.  until  the  armature  takes  her 
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‘bite’  from  the  rail.  In  interpreting  Fig.  2,  it  should  be  remembered  that  the 
waves  shown  are  standing  waves,  stationary  with  respect  to  the  armature.  This 
process  is  repeated  for  the  rest  of  the  passage  of  the  armature  along  the  gun. 
What  determines  the  gouging  frequency  is  not  yet  explained,  but  as  the  armature 
continues  to  accelerate  and  9  increases,  the  wave  length  decreases,  so  we  expect 
the  gouging  frequency  to  increase.  Also  the  amplitude  of  the  waves  decreases, 
so  if  the  model  is  correct  we  should  see  a  reduction  in  the  depth  of  the  gouges 
along  the  rail.  to  b  ^ 


3.  Response  of  the  foundation 


The  most  difficult  parameter  to  determine  in  the  above  theory  is  the  amplitude 
of  the  spring  constant  F.  In  fact  the  foundation  can  only  be  roughly  approxi¬ 
mated  by  a  series  of  independent,  local  springs.  Our  aim  here  is  to  find  a  good 
approximation  to  F. 


armature 


copper  rail 
(beam)  \ 


L  beam  load 

V 


spread  of  load  onto 
foundation 
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i 

1 

b 


Figure  3:  A  rail  gun  beam  and  foundation 

In  Fig.  3  we  give  an  overall  picture  of  the  load,  beam  and  foundation,  with 
the  important  dimensions  identified.  It  is  considerably  out  of  scale.  Figure  4 
shows  the  compression  of  the  foundation  by  a  uniform  load  L  applied  to  the 
beam  over  an  infinite  strip  of  width  a.  This  load  is  transferred  through  the  beam 
to  the  foundation.  One  problem  that  needs  consideration  is  the  distribution  of 
pressure  between  the  beam  and  the  foundation.  To  avoid  solving  a  complicated 
boundary  value  problem,  at  this  stage  we  shall  make  a  rough  approximation  by 
assuming  that  the  load  is  equivalent  in  its  effect  to  a  uniform  pressure  p  acting 
on  the  foundation  over  a  range  2 h,  where  the  relation  between  the  length  a  of  the 
armature  and  the  distance  h  is  to  be  determined. 
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For  h  we  shall  adopt  the  approximate  relation 

2h  =  a  +  b.  (15) 

This  means  that  the  spread  on  each  side  of  the  load  is  equal  to  half  of  th* 

the^ext  bv  T°  ™S  empiricaI  relati°n  follows  from  Fig  34  in 

the  text  by  Timoshenko  and  Goodier  [31  which  fimIrP  t  *  w,  r  S' , 

paper  written  in  1903  by  Filon  [4].1  ’  r0m  a  remarkable 


r  igure  4:  btrain  of  the  foundation 

rigid  base  on  y  =  0  and  a  uniform  load  onvL  1  l  S  Jt  has  a 

beam,  over  a  distance  2 h.  The  load  is  directed  donft^ga^^O^  °n  W 
can  replace  the  rigid  boundary  by  mirroring:  the  load  in  a  •  m  *  7  °Y  We 
the  beam,  for  evidently  the  axis  of  such  a  System  is  strakte  fl  bT^  ben“fth 
were  rigid.  We  require  the  displacement  ,  of  the  surSe  u  =  Welhlh  ?? 
this  by  integrating  the  stress  strain  relation  h  deduce 

cfo  _  l 

dy  ~  Yf  ~  u^x  +  '*))  -  (16) 


We  introduce  the  dimensionless  parameters 


to  a  distribution  of  forced  olfa^?  sides ' blSe^of1hemeaun?  the  StreSS  “  *  parallelopiPed’  subject 
their  presentation.  0me  °f  the  equat,ons  req^ed  more  than  one  page  for 
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deflection  v 


Figure  5:  Laminar  subject  to  symmetrical,  concentrated  loads 
then  find  that  the  stresses  are  given  by  (cf.  pp.  48-50  of  [3]) 


ay  = 


—  _L  j  cosh  ma  +  sinh  ma 

l  ^  m= i  |  sinh  2 ma  4-  2mo:  c°sh  (may/c) 


(may/c)  sinh  (may/c)  sinh  ma ) 

sinh  2ma  +  2ma  /  cos(m^/h) .  (18) 

ai  is  necessary  only  to  change  the  signs  of  the  indicated  symbols 
Integrating  over  0  <  y  <  c,  we  find 


pc  h  pc 

v  ~eJ  , 

where  $(X:  a,  3)  =  -  "V  -1  (  sinh2  sin  Pm  \  ,  „x, 

*  m=i  am2  \sinh  2 ma  +  2ma )  cos(m^)- 


(19) 


We  find  from  (19)  that  v  changes  sign,  as  indicated  in  Fig.  5  at  x  -  ±kh 
negative  ^  4  “*  5  ^  -It 

H  =  TkLvd{x^h^ 

We  shall  therefore  define  F  by 

F\v\  =  L/2h,  L  =  2hwrp ,  ^20) 

wtae  u,  is  the  width  of  the  rail  (see  Fig.  3).  Since  we  are  ignoring  stress  in 
the  OZ- direction,  it  does  not  matter  whether  the  foundation  is  regarded  as  beine 
a  laminar,  or  an  infinite  collection  of  strips  of  width  wr,  comprising  the  lamf 
nar  The  quantity  F  is  the  restoring  force  acting  on  the  beam  due  ton  uniform 
displacement  v  of  the  foundation.  This  allows  us  to  regard  F  a«  k  •  ,  ? 

‘spring’  force.  It  follows  from  (19)  and  (20)  that  §  b  S  *  ^ ’ 

jpi _  wrE 

~  c{(h/ /)  +  ($>}  ’  (21) 
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where 


<4>  =  S  /*♦(*/*; 


..  yhe,re  “  °ne  further  Wtoximation  that  should  be  mentioned.  We  have  placed 
the  load  at  the  center  of  the  beam  to  simplify  the  calculation;  of  course  gouging 
w.11  not  necessarily  occur  at  this  position,  but  because  h/l  is  typical  ve™  small 
one  expects  that  the  resulting  error  will  be  negligible.  ’ 


4*  Application  of  the  theory 

The  following  values  are  taken  from  the  rail-gun  experiments  at  1AT,  University 
of  Texas  at  Austin.  The  dimensions  shown  in  Fig.  3  are:  7 

a  =  18.3  mm,  6  =  6.35  mm,  c  =  50.4  mm,  wr  =  31.75  mm,  /  =  600mm 

!°  - ^  haV6:  Ef  =  1010  Pa>  =  1.68x10“  Pa 
Pf  -  1900 kgm  ,  pcu  =  8960kgm-3,  pal  =  2700 kgm"3,  p  =  1.81  kgnr2. 

T  tFe_S,e0Va4IUeS  and  ^  and  (17)  we  §et  a  =  0-266,  0  =  0  0647 
a  7JcZ  °ofXPh.  ’  ^  =  H4ntm2  ^  f0U°Wing  table  SetS  «*  values  of  *  as 

Table  1:  The  function  <&(x/h) 


x/h 

0 

0.2 

0.4 

0.6 

0.8 

1.0 

1.2 

1.4 

1  6 

]  R 

103$ 

621 

615 

597 

565 

513 

416 

316 

257 

212 

x  .0 

176 

x/h 

2.0 

2.2 

2.4 

2.6 

2.8 

3.0 

3.5 

4 

4.5 

5  n 

103$ 

146 

121 

99 

81 

65 

51 

24 

5 

-7 

U  .\J 

-16 

We  find  $  0  at  x/h  «  4.25  and  that  its  average  over  (0  4  25  h)  is  n  oon  tv. 

from  (21)  we  get  F  =  2J110-.  Finally  the  cola,  spe^loUotl  from  (1)  Jt? 

v.-fi™)1. 


[  P  7 


(22) 


detemintag'L  Iponst’  of' found”  tWst  Xftotte"8  inVOlV<i<i  “! 
"o^u^  -  P^y^'a^ole. 

against  observation,  before  the  theo^  Z  — 
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5.  Copper  armatures 


The  above  theory  appears  to  have  some  merit,  since  it  offers  an  explanation 
of  why  gouging  might  occur  and — for  Aluminium  armatures — arrives  at  a  value 
of  the  critical  speed  in  close  agreement  with  observation.  There  remains  one 
problem  to  consider.  When  the  armature  is  made  of  Copper  and  all  the  other 
conditions  remain  the  same,  gouging  sets  in  at  speeds  of  only  600  m/s,  well  below 
the  value  obtained  above.  There  are  two  possibilities:  either  the  theory  is  wrong, 
or  that  Copper  armatures  behave  differently,  perhaps  by  exciting  antisymmetrical 
displacements,  as  indicated  in  Fig.  6. 


rails  moving  past  armature 


Figure  6:  Antisymmetric  displacements  of  the  rails 

Because  its  mass  per  unit  length  is  greater  than  than  of  the  rails  supporting 
it,  a  Copper  armature  is  vulnerable  to  small  modulations  in  the  straightness  of 
the  rails.  It  is  therefore  less  likely  to  travel  smoothly  along  the  centerline  of  a 
rail  gun  than  is  a  much  lighter  Aluminium  armature.  Transverse  oscillations  of 
the  type  indicated  in  Fig.  6  may  be  excited,  and  if  the  armature  speed  reaches 
a  critical  value,  gouging  becomes  possible.  We  shall  accept  this  description  and 
determine  a  value  for  Vc,  which  we  shall  identify  with  V*. 

We  now  regard  the  armature,  labelled  A  in  Fig.  6.  as  being  an  additional 
concentrated  mass  M,  on  the  beam  and  apply  Newton;s  second  law  to  it.  As  A 
moves  between  the  rails,  first  one  rail  then  the  other  provides  the  restoring  force. 
There  are  two  forces  acting  on  A,  namely  the  difference  in  shear  across  A  in  the 
rail  that  is  decelerating  the  mass  and  the  force  due  to  the  foundation  behind  that 
rail.  The  equation  of  motion  for  A  is  therefore 


-El 


(23) 


where  [•  •  • }A  denotes  the  jump  in  •  •  •  and  p  is  the  mass  of  the  armature  plus  rails 
per  unit  length.  In  order  to  obtain  a  rough  estimate,  we  shall  write 


i[£vl  „(#v\ 

a  dx*  ~  Urr4 
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which  is  exact  in  the  limit  as  a  -*•  0.  Now  (23)  assumes  the  form 

r -,r93 4y  „<92y 

EId?  +  pW+Fy  =  0’  (24) 

which  is  identical  with  (5),  with  the  modifications  that  p  replaces  p  and  that  the 
equation  applies  to  only  to  the  armature  and  not  to  the  whole  rail-gun  system. 
Despite  this,  it  should  give  an  estimate  of  the  critical  velocity  for  antisymmetrical 
displacements. 

The  theory  is  now  the  same  as  in  §1,  leading  to  equation  (10): 


4  FEI\* 

P2  )  ’ 


(25) 


where  p  =  2p+M/a. 

The  Aluminium  armatures  in  the  experiments  cited  in  the  previous  section  have 
a  mass  of  35.31  grm.  A  similar  Copper  armature  would  therefore  have  a  mass  of 
M  =  117  grm.  Hence,  with  the  values  qiven  in  §4,  p  =  2p  +  M/a  =  10.0  kg  m-2 

It  follows  from  (25)  that  V*  =  587 ms-1,  in  close  agreement  with  Dr.  Marshall’s 
observations. 

A  simple  test  of  this  theory  would  be  to  load  an  Aluminium  armature  with 
heavy  metal  and  see  if  this  reproduces  the  value  of  V*  for  a  Copper  armature 
And  other  tests  are  suggested  by  the  functional  dependence  exhibited  in  (22)  and 
(25).  If  these  tests  confirm  the  theory,  then  the  more  difficult  task  of  modeling 
the  gouging  process  itself  would  the  next  step. 
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Three-Dimensional  Effects  on  the  Melt-Wave  Model 

Summary 


During  my  remit,  visit  (26th  April  1 7 1. 1 1  May,  1996)  to  IAT,  I  was  shown 
a  recovered  armature,  with  a  marked  three-dimensional  pattern  of  erosion,  as 
indicated  in  Fig.  1  below  the  contact  surface  of  the  armature  was  deeply  wasted 
along  its  edges.  To  interpret,  this  as  due  to  a  melt  wave  advancing  from  the  rear 
of  the  armature  would  require  a  model  in  which  the  speed  of  the  wave,  vm.  was  a 
function  of  the  distance  2,  measured  orthogonal  to  the  direction  of  propagation. 

The  speed  vm  is  strongly  dependent  on  the  total  current  I  flowing  across  the 
interface  between  the  armature  and  the  rail.  In  the  two-dimensional  model,  this 
is  treated  as  being  independent  of  2.  Let  i]r  be  the  rail  resistivity,  V  the  speed 
of  the  armature  and  //,,  the  permeability  of  free  space,  then  fir  =  rjr/(/i(]V)  is  the 
very  small  diffusion  length  of  the  current,  measured  in  the  direction  of  motion. 
If  w  denotes  the  depth  of  the  surface  of  the  armature  in  contact  with  a  rail,  the 
total  current  is  I  =  Srwj0 ,  where  jo  is  the  current  density  immediately  in  front  of 
the  melt-wave.  The  method  adopted  below  to  introduce  three-dimensional  effects 
into  the  melt- wave  model,  is  to  replace  I  by  j0  in  the  expression  for  vm ,  and  then 
to  obtain  an  estimate  of  the  distribution  jo{z),  using  approximate  diffusivity 
arguments.  This  approach  allows  most  of  the  previous  melt  wave  theory  to  be 
retained. 

The  report  begins  by  describing  a  small  numerical  correction  to  the  two- 
dimensional  model  presented  in  [1],  Then  an  expression  for  the  minimum  current 
required  to  initiate  the  melt  wave  and  hence  to  provide  the  lubrication  to  allow 
the  armature  to  move,  is  derived. 

To  improve  on  the  diffusivity  approach,  it  is  necessary  to  undertake  a  four¬ 
dimensional  (space  and  time)  numerical  calculation.  A  method  of  doing  this, 
that  does  not  involved  more  than  three  independent  variables  at  each  step,  is 
described  at  the  end  of  this  report. 


Figure  1:  Erosion  of  armature  due  to  three-dimensional  currents 
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1.  Some  remarks  about  the  two-dimensional  model 


In  a  correction  made  to  the  version  of  [1]  presented  at  Baltimore,  I  have  changed 
the  value  of  the  constant  n  appearing  in  the  paper  from  an  incorrect  value  of 
1.148  to  unity.  The  paragraph  including  equation  (9)  of  that  paper  should  now 
he  replaced  by: 


“....the  sum  over  the  separation  parameter  rn  of  terms  of  the  type 

:i'y  =  Am  cos (my/6r)eax/6r ,  (9) 

where  «  =  ^{(1  +  4m') 2  +  1}  (3;  <  0,  y  <  0). 


To  determine  rn  it  is  necessary  to  consider  the  diffusion  of  the  electric 
current  into  the  rail  as  viewed  from  the  armature.  In  a  time  interval  6t ,  a 
current  element  at  (xo,0)  diffuses  a  distance  y  =  —  into  the  rail  as 

it  is  convected  a  distance  x  =  VSt  along  OX.  The  element  therefore  moves 
along  the  parabola  y 2  =  6r(x  —  no).  Thus,  if  no  =  —fc2a,  0  <  k  <  1,  is  the 
distance  between  the  melt  wave  and  the  front  of  the  armature,  y  =  k(6ra)? 
is  the  value  of  y  at  which  jy{0,  y)  is  zero.  Hence,  by  (9),  the  smallest  value 
of  m  is  mo  —  ir(6r/a)/2k.  This  corresponds  to  the  smallest  value  of  a, 
and  hence,  for  a  given  (negative)  value  of  n,  it  gives  the  dominant  term 
in  jy(x,  0).  As  6r/a  is  typically  1CT3,  except  for  very  small  values  of  k , 
mo  ~  0  and  a  ^  1,  which  value  we  shall  adopt.” 

With  the  modification  just  described,  the  velocity  vm  of  the  melt  wave  and 
the  depth  h  of  the  material  melted  (and  hence  presumably  removed  by  adherence 
to  the  rail),  is  given  by 


VJIqPV  2jxqGV 
2rjrpQw 2  nr]r 


and 


TTTfr 

2p0V' 


where 


Q  =  c; >{T,n  ~  To)  +  Lm  , 


(3) 


In  these  expressions  Tm  is  the  melting  temperature,  T()  is  the  laboratory  temper¬ 
ature,  7)a  is  resistivity  in  the  solid  armature,  averaged  over  ( Tm ,  T0),  r/r  is  the  rail 
resistivity,  /  is  the  total  current  flowing  through  the  armature,  //.0  is  the  perme¬ 
ability  of  free  space,  p  is  the  armature  density,  w  is  the  depth  of  the  surface  of 
the  armature  in  contact  with  a  rail,  Lm  is  the  enthalpy  of  fusion  (latent  heat), 


46 


<)>  is  l lie  specific  heat  at  const ;u 1 1  pressure.  \„,  \r  are  the  thermal  diffusivities  in 
the  armature  and  rail,  and  \  is  the  speed  of  the  armature. 

The  total  current  I  in  (2)  is  related  to  the  current  density  into  the  armature 
at  the  interface  between  the  armature  and  rail,  sav  j(),  bv 

I  =  »'<y/o  ■  (4) 


Of  the  variables  in  (2).  the  only  one  that,  could  depend  strongly  on  2  is  the 
current  density.  j0  implicit  in  /.  It  follows  that,  in  order  to  model  the  erosion 
shown  in  Fig.  1,  we  must  determine  the  2  dependence  of  the  current  density  across 
the  interface. 


Notice  from  (2)  that  the  melt  wave  will  not  begin  until  a  critical  current  given 
by 


/*  =  Kw , 


(5) 


where 


4  pQG'V 

*Va  ) 


is  attained.  For  copper  rails  and  aluminium  armatures,  we  find  the  values  (in 
SI  units):  G  =  0.589,  fja  =  7.3xlCT8,  *0  =  1.07xl(T4,  Xr  =  1.16><10_4,  p  = 
2.66xl03,  Q  =  9.67xl05.  Substituting  into  (4),  we  get  K  =  2.37xl06.  For 
example,  with  w  =  18.85xl0_3m,  the  critical  starting  current  is  I*  =  4.48xl04  A. 
If  it  is  necessary  to  lubricate  the  rails  with  liquid  aluminium  before  it  will  move, 
then  (5)  gives  the  current  necessary  to  start  the  armature  in  motion.  It  would  of 
interest  to  use  IAT  experiments  to  check  this  speculation. 


2.  An  semi-empirical  treatment  of  the  3-D  current 

Figure  2  below  shows  a  cross-section  of  a  rail  just  behind  an  armature.  The  shaded 
area  of  this  section  indicates  the  region  through  which  the  current  is  flowing;  the 
closer  it  is  to  the  armature,  the  more  this  region  will  be  distorted  upwards.  Our 
task  is  to  obtain  an  estimate  of  the  current  distribution,  jy(x,0,  2),  through  the 
interface  on  y  =  0,  between  the  rail  and  the  armature.  In  the  two-dimensional 
description  of  the  melt  wave  given  in  [1],  we  had 

jy(x,  0)  =  jl)ex/6r  (x  negative) .  (6) 

The  only  accurate  method  of  determining  a  3-D  generalization  of  (6)  is  via  a 
numerical  solution,  which  requires  time  to  be  included  as  an  independent  variable. 
We  shall  outline  a  method  of  doing  this  in  the  following  section.  We  should  try 
to  avoid  this  complication  for  each  application,  although  some  calculations  will 
bo  necessary  to  support—  or  perhaps  to  eliminate-  -the  empirical  approach  given 
here.  The  main  difficulty  in  finding  an  analytical  extension  of  (6)  is  that  both 
flow  and  diffusion  contribute  to  the  three-dimensional  effects. 
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Figure  2:  Three-dimensional  currents  across  the  rail /armature  interface 

Our  method  is  to  allow  for  diffusion  by  a  factor  (1  +  (3 exp[|z|/7]).  where  the 
parameters  jd  and  7  are  to  be  determined.  Thus  our  hypothesis  is  that  the  normal 
current  through  the  interface  is  given  by  an  expression  of  the  form 

jy(x,0,z)  =  cr(l  +  (iSz^'r)ec/Sr .  (7) 

The  total  current  I  passing  through  the  armature  is  given  by 

rwj 2  rO 

I  =  2  dz  jy  dx , 

JO  7— 00 

where  w  is  the  width  of  the  armature  (see  Fig.  2).  From  (7) 

/  =  Sra[w  +  2/Fy(e’"/27  —  1)] .  (8) 

In  Fig.  3  we  have  sketched  the  current  distribution  at  a  fixed  value  of  x.  The 
2-D  distribution  is  flat,  with  the  same  area  beneath  it  as  the  3-D  case.  The  skin 
depth  in  a  section  immediately  below  the  armature,  and  at  the  position  of  the 
melt  wave,  is  given  by  2A,1  where 

A  =  (&•*)*  =  (W'O*  (fr  =  Vr/lMi)  ,  (9) 

in  which  £r  is  the  magnetic  diffusivitv  in  the  rail  and  a  is  the  length  of  the 
armature  in  contact  with  the  rail.  The  time  t  is  the  armature  transit  time,  since 
this  is  the  time  during  whivli  the  rail  is  exposed  to  the  current.  The  factor  ‘2’ 
takes  care  of  the  flow  of  the  current  layer  from  beneath  the  rail,  as  indicated  in 
the  figure. 

‘In  copper,  at,  room  temperature,  A  «  O.llG/i,  while  at  the  melting  temperature  of  alu¬ 
minium,  A  «  0.219/3 . 
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Figure  3:  The  current  density  distribution 


We  now  assume  that  the  current  flowing  up  through  the  rail  at  a  given  point 
z  is  proportional  to  the  depth  of  the  rail  carrying  current,  just  before  the  flow 
turns  “upwards”  into  the  armature.  Let  js,  jc  be  the  surface  and  central  current 
densities,  as  shown  in  Fig.  3,  then 

Js  f  b  /in\ 

~  =  J  7tT  =  9 1  say,  (10) 

Jc 

where  b  is  the  thickness  of  the  rail  and  /  is  a  form  factor  to  allow  for  the  dis¬ 
tribution  of  the  skin  depth  on  the  side  of  the  rail.  We  expect  this  to  be  roughly 
triangular,  with  a  thickness  of  A  at  the  top.  This  suggests  a  value  f  =  \- 

It  follows  from  (7)  and  (10)  that 


P  = 


P-1 

exp{w/2j)  -  g 


(1  <  g  <  ew,2y). 


(11) 


There  remains  the  parameter  7.  We  shall  assume  that  the  edge  current  diffuses 
inwards  and  set  7  equal  to  A. 

Summarizing,  we  now  have 


jy(x,  0,  z)  =  o:(l  +  (Je^/x)ex/6r  =  jo€x/Sr , 

/ 


where  a  = 


7  /?  = 


P-1 


(12) 


8t{w  +  2/lA[exp(u>/2A)  —  1]}  ’  exp(w/2A)  —  1 


To  apply  the  above  theory  to  a  railgun,  we  omit  the  second  right-hand  term 
in  (2),  since  it  is  much  smaller  than  the  first.  Then  from  (2).  and  (4)  we  get 


,2 

4/;Qh  70 
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which  is  independent  oi  the  width  v-  oi  the'  annaturr.  It  is  therefore  an  expression 
we  can  apply  locally.  It.  follows  from  (12)  that 

=  T7^',2(1+rfe':l/A)2-  (13) 

\/)Qh  v  ’ 

The  2-D  result  is  recovered  at  d  =  0. 

The  dependence  of  vm  on  2  means  that  the  melt-wave  front  will  develop  edge 
spurs,  as  shown  in  Fig.  1.  The  reason  why  we  do  not  assume  that  the  melt  wave 
propagates  inward  from  the  spurs  is  that  the  length  scales  A  and  br  appearing  in 
(12)  have  very  different  magnitudes,  which  means  that  the  current  is  concentrated 
at  the  smallest  values  of  .7;,  so  the  melt  wave  proceeds  in  the  OX-direction  rather 
than  in  the  OZ-direction. 

What  is  now  required  is  a  numerical  evaluation  of  jy(x,0,z)  in  typical  condi¬ 
tions,  in  order  to  check  the  structure  of  the  first  equation  in  (12).  In  the  following 
section,  we  shall  outline  a  numerical  method  of  doing  this. 

Illustration  (using  I  AT  values) 

The  temperature  of  the  rail  in  the  vicinity  of  the  armature  is  required  for 
an  accurate  determination  of  A.  We  shall  take  the  average  of  the  values  at  room 
temperature  and  aluminium  melting  temperature,  quoted  in  the  footnote  on  page 
5.  At  a  =  18.3mm,  and  V  =  500 ms-1,  (9)  gives  A  =  1mm.  Take  /  =  |  in  (10). 
With  b  —  6.35mm  (|w),  we  get  g  =  1.59,  i.e.  js/jc  =  1-59  (see  Fig.  3).  And 
with  w  =  18.85  mm,  (12)  gives  (3  =  0.59  exp(— w/2X)  =  4.71xl0~5.  The  current 
distribution  is  therefore 

jy(0,0,2)  =  cv(l  +  0.59exp{(|z|  -  |io)})  (distances  in  mm). 

The  melt-wave  velocity  is 

vm  oc  (1  +  0.59 exp{(|z|  -  \ w )})2  . 

The  speed  on  the  edge  of  the  armature  is  1.592  =  2.53  times  faster  than  in  the 
centre.  Of  course,  these  values  are  merely  illustrative,  since  the  factor  /  in  (10) 
is  really  unknown  and  the  ratio  just  calculated  depends  om  /2. 

3.  Numerical  treatment  of  the  3-D  current 


In  Fig.  4  we  show  a  simple  example  of  a  railgun — one  rail  and  half  an  armature, 
so  that  E  is  a  surface  of  symmetry.  The  armature  is  at  a  distance  s  =  Vt,  where 
t  is  the  elapsed  time,  more  accurately, 

,s  =  I'  V  dt ,  ds  =  V  dt. .  (14) 

./o 

It  is  convenient,  to  use  a  in  place  of  tin;  time  t. 
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Figure  4:  One  rail  and  half  the  armature 


The  breech  takes  the  initial  current  through  a  thin  skin  depth.  For  a  small 
value  of  s,  say  so  corresponding  to  t  =  f0>  this  is 

Ao  =  (&«<>)*.  (15) 

At  t  =  0,  the  armature  is  stationary  at  the  breech  and  the  current  density  is 
momentarily  infinite.  It  is  therefore  necessary  to  delay  slightly  before  starting 
the  calculation.  At  t  =  to,  the  cross-sectional  area  available  to  the  current  is  (see 
Fig.  4)  2As(tw  +  h),  so  the  initial  current  density  has  the  components 

jx(to)  =  -f(^o)/ {2A0(u;  +  h)},  jy(to )  =  0,  jz{to)  =  0 .  (16) 


The  boundary  conditions  for  j  are  evident — no  currents  normal  to  the  rail 
surface  and  symmetry  across  the  surface  E. 

Relative  to  a  frame  fixed  in  the  armature,  the  rail  has  a  velocity  along  OX  of 
V,  for  which  case  the  curl  of  Ohm’s  law  yields 

dB  <9B  , 

~dt+V!h~^V  ' 

Since  /ioj  =  VxB,  the  curl  of  this  equation  yields 


m  =  v^  + 


By  (14)  we  have 


dj  ,  0j 


+ 


dx  ds  ’ 


dx  dt 


(6  s  {,/V  =  «(»)) . 


Given  V(t),  we  find  s(t)  =  /0£  V(t)  dt,  and  invert  this  to  get  t(s),  then  obtain 
V(t(s)).  We  must  also  ensure  that 

V-j  =  0.  (18) 
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The  domain  shown  in  Fig.  1  may  be  halved  bv  using  the  symmetry  about  the 
O  A” -axis.  Also  it  may  be  .sufficiently  accurate  to  terminate  t  he  domain  about  two 
armature  lengths  beyond  the  leading  edge  of  the  armature. 

The  armature/rail  interface  is  the  surface  where  (see  figure),  T,  =  (—a  < 
x  <  0,  ij  =  0.  —r,u)  +  (■  <  z  <  r,ii!  —  ().  The  task  is  to  determine  y,  z.  a). 
in  particular  to  find  |y|  on  S,.  The’  boundary  conditions  should  ensure  that  the 
total  current  entering  the  breech,  also  crosses  S:  this  requirement  will  be  a  useful 
check  on  the  calculations. 

Four  independent  variables  are  one  too  many  for  practical  calculation.  The 
time  variable  s  may  be  accommodated  by  the  following  iterative  scheme.  Start 
with  the  initial  steady  state, 

<5°v2j°  =  s  =  .so,  V'  =  Vo,  So  =  tr/Vo , 

then  advance  in  small  steps  of  s: 

c  V72  *  1  |  J 1  Jo  i  A 

«3iV  ji  —  H — — — ,  s  —  so  +  As, 

. . .  = 

c  v72i  ^Jti  Jn  Jn — 1  .  a 

SnV  jn  =  H - 2^ — ’  s  =  s0+nAs, 

continuing  until  s  equals  the  length  of  the  rails.  During  each  step,  the  condition 

djxn  fyjyn  djzn  ^ 

dx  dy  ~dz'~  ’ 

must  be  satisfied. 


Reference 


1.  Woods,  L.C.  “The  current  melt  -wave  model”.  IEEE  Trans,  on  Magnetics , 
Vol  32,  No.  1, .  1997. 


Progress  Peport  on  Research  into  Electromagnetic  Launchers 


Boundary  Conditions  Across  a 
Pail  /  Armature  Contact  (Surface 


June  1996 


Professor  L.  C.  Woods 
University  of  Oxford 
Mathematical  Institute 


Boundary  Conditions  Across  a  Rail  /  Armature  Contact  Surface 

Summary 

In  this  report  I  shall  give  an  account  of  the  conclusions  that  I  arrived  at  during 
my  visit  to  IAT  in  May,  concerning  the  averaged  boundary  conditions  appropriate 
for  a  contact  surface  comprised  of  “a-spots” .  One  aim  was  to  obtain  a  macroscopic 
value  for  the  resistivity,  say  rjeff,  that  could  be  used  across  the  interface  between 
the  rail  and  armature  in  place  of  the  “microscopic”  details  of  the  a-spots.  This 
was  followed  by  a  calculation  of  the  heat  generated  in  the  contact.  The  effect  of 
relative  motion  between  the  rail  and  armature  was  then  estimated  and  applied 
to  the  melt- wave  model  [1]. 

Since  my  visit,  I  have  reconsidered  the  a-spot  model  and  concluded  that  in  the 
presence  of  a  high  speed  of  relative  motion  between  the  surfaces,  it  has  little  direct 
relevance  to  the  average  boundary  conditions  across  the  interface,  except  perhaps 
at  the  very  beginning  of  the  armature  motion.  I  have  therefore  developed  an 
‘asperity-plough’  model,  in  which  copper  asperities  plough  through  the  armature 
material,  smearing  contact  points  that  initially  behave  momentarily  like  a-spots, 
into  narrow  streaks,  stretching  from  the  site  of  the  initial  spot,  backwards  to  the 
rear  of  the  armature. 

Provided  the  surfaces  remain  pressed  together,  the  intense  local  heating  at  the 
asperities  melts  sufficient  aluminium  to  allow  the  streaks  to  run  together  to  form 
a  continuous  electrical  contact  between  the  rail  and  the  armature.  Thus  we  arrive 
at  a  model  very  similar  to  the  melt-wave  model,  except  that  the  ohmic  heating 
is  more  intense  due  to  the  initial  current  concentration  and  the  depth  of  material 
entrained  on  the  rail  depends  on  the  size  of  the  asperities,  or  equivalently,  on  the 
thickness  of  the  contact  layer  and  not  on  the  skin  depth  of  the  electric  current. 

At  the  end  of  the  report,  we  offer  an  explanation  of  the  transverse  stria- 
tions  visible  on  the  entrained  material,  when  both  the  rails  and  the  armature  are 
made  of  aluminium.  It  is  very  likely  to  be  due  to  a  shear-layer  instability,  such 
as  that  observed  in  the  wing-tip  streamers  generated  by  high-flying  aircraft  in 
supercooled,  atmospheric  conditions. 


1.  The  asperity-plough  model 


There  are  two  distinct  types  of  asperity — those  consisting  of  copper  intrusions  into 
the  aluminium  armature  (c-asperities)  and  the  inverse  type,  i.e.  those  formed  by 
aluminium  intiusions  into  the  copper  rails  (a-asperities).  Figure  1  illustrates  a 
c-asperity.  We  may  suppose  that  initially,  both  types  are  equally  numerous— of 
course  this  will  depend  on  the  smoothness  to  which  the  surfaces  are  finished. 
We  would  expect  a-asperities  to  behave  rather  like  a-spots;  they  move  with  the 
armatuxe  and  aie  feed  with  current  from  the  line  of  points  immediatelv  opposite 
on  the  rail.  However,  there  are  two  reasons  why  the  lifetime  of  a-asperities  must 
be  much  shorter  than  that  of  c-asperities.  First,  at  high  temperatures,  aluminium 
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has  only  a  fraction  of  the  shear  strength  of  copper  and  secondly  it  melts  at  a  much 
lower  temperatures  than  copper.  Thus,  relative  to  c-asperities,  a-asperities  are 
sheared  off  quickly  and  even  if  they  survive  until  melting,  at  this  stage  they 
immediately  disappear  by  being  entrained  on  the  rail.  This  description  implies 
that  the  current  into  the  armature  is  largely  transferred  through  c-asperities,  or 
rather  through  the  strips  of  molten  aluminium  generated  by  these  asperities,  as 
they  plough  through  the  armature  during  their  transit  towards  its  trailing  edge. 
The  distinction  we  are  advancing  between  the  two  kinds  of  asperity  receives  some 
support  by  the  observation  that  with  aluminium  armatures  on  aluminium  rails, 
the  appearance  of  the  material  entrained  on  the  rails,  with  strong  transverse 
striations,  is  very  different  from  the  standard  copper  rails  (see  §5  below). 

armature 


Figure  1:  A  copper  asperity  at  an  early  stage 

In  Fig.  1  we  show  a  copper  asperity,  C,  at  the  stage  when  the  current  it  carries 
has  melted  the  surrounding  aluminium.  Because  of  the  velocity  skin  effect,  the 
dominant  current-carrying  asperities  will  begin  their  lives  close  to  the  trailing 
edge,  TE,  of  the  armature.  Initially  the  current  will  be  concentrated  near  the 
surface  between  C  and  TE,  so  there  will  be  some  melting  in  front  of  C.  As  C 
moves,  it  will  leave  behind  it  melted  aluminium,  but  this  will  soon  freeze  on  to 
the  rail,  leaving  a  high-resistance  gap  trailing  it.  Therefore,  during  the  transit  of 
C  from  its  initial  site  to  the  trailing  edge,  there  will  be  good  electrical  contact 
across  a  narrow  strip  extending,  either  from  C’s  initial  site  to  TE,  or  from  a 
freeze-wave  just  behind  C  to  TE.  When  C  passes  beyond  TE,  the  strip  will  be 
entrained  on  the  rail  and  carried  away.  This  is  illustrated  in  Fig.  2.  However,  it 
should  be  remarked  that  the  figure  is  not  even  close  to  scale,  since  the  spots  are 
believed  to  be  typically  about  1/10  mm  in  diameter,  with  TE  some  two  hundred 
times  larger. 

As  the  armature  accelerates  along  the  rails,  more  and  more  c-asperities  are 
'activated’,  until  the  strips  they  generate  form  a  near-continuous  region  of  good 
conductivity  near  the  trailing  edge  of  the  armature.  Provided  the  armature  and 
rails  continued  to  be  pressed  together,  new,  albeit  shorter,  asperities  will  be 
continue  the  process.  But  as  the  material  near  the  rear  of  the  armature  is  lost 
to  the  rails,  the  effective  trailing  edge  will  move  forwards,  in  much  the  same 
way  as  the  melt  wave  in  the  ideal  contact  model.  On  this  model,  we  expect  the 
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Figure  2:  Conducting  streaks  generated  by  c-asperities 

aluminium  deposited  on  the  rails  to  have  a  striated  appearance;  this  can  be  seen 
in  photographs  of  the  rails  [2]  (also  see  Fig.  8). 


2.  Some  results  for  the  a-spot  model 

Since  in  the  initial  stages,  the  contact  between  the  armature  and  rail  will  be 
through  a-spots,  we  shall  start  by  reviewing  the  theory  of  these  spots.  In  the 
following,  we  shall  adopt  the  subscript  ‘s’  to  indicated  values  for  a  typical  a-spot; 
quantities  without  a  subscript  will  apply  to  the  ideal  (perfect  contact)  values. 

armature 


... 

j 

c-asperities  — 

Figure  3:  Idealized  interface  of  c-asperities 

Let  As  be  the  total  area  of  solid  or  liquid  contact  between  the  armature  and 
tail  and  A  the  actual  or  nominal  area  of  the  contact  surface,  then  the  ratio 
^  —  As/u 4  will  be  much  less  than  unity — values  of  0.1  to  0.2  are  suggested  in  the 
literature.  Let  F  denote  the  normal  pressure  between  the  rail  and  armature,  then 
the  nominal  pressure  is  p  =  F/A.  For  given  metals,  e  will  depend  on  p  and  also 
on  the  surface  temperature.  Thus 

e  =  As/A,  e  =  e(p,  T) .  (1) 
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If  I  is  the  total  current 


js  =  I /As,  j  =  I /A,  so  js  —  j/e. 


(2) 


Let  Tjs  be  the  spot  resistivity,  then  (see  Fig.  3) 

l  —  a  a 
Vs  =  —j—Va  +  y  Vr  , 


(3) 


where  a  is  the  average  ‘radius’  of  the  asperities,  /  is  the  thickness  of  the  interface 
layer,  and  rja  and  rjr  are  the  resistivities  of  the  armature  and  rail.  We  shall 
assume  that  only  the  armature  material  is  melted.  At  the  melting  temperature 
of  aluminium  (933  K),  rja  =  10.7xl0~8nm  and  rjr  =  6.3*10-8Qm.  The  ratio 
a/l  is  not  easy  to  estimate.  Since  r]s  lies  between  r]r  and  7]a  and  these  values 
are  not  wildly  different,  we  shall  be  content  to  accept  the  average  spot  value, 
Vs  =  8.5xl0-8  m. 

Let  E  denote  the  electric  field  across  the  layer,  then 


Vsjs  =  E  —  AV/l ,  AT f  =  IR,  R  =  Vsl/As, 


where  R  is  both  the  layer  resistance  and  the  spot  resistance.  The  effective  resis¬ 
tivity  for  the  layer  therefore  follows  from 


p  _  Vsl 


^7eff  5 


Vex  =  Vs/e . 


(4) 


Let  Ns  be  the  spot  number  density,  then 


N.  = 


A„ 


it  a2  A 


—  •  L  =  N-1' 2 

_o  !  s  ) 


ITQ,* 


where  a  is  the  average  spot  radius  and  Ls  is  the  spot  spacing.  Thus 

9 

7 ra 


L2 


(5) 


(6) 


To  illustrate  the  theory,  we  shall  adopt  [3]  a  =  0.5xl0-4  m  and  Ls  =  2.3xl0~4  m, 
which  gives  e  «  0.15  and  Ns  «  2xl07m-2.  Hence  from  (4)  and  the  estimate  for  r]s 
given  above,  7?eff  =  5.7xl0_7flm.  The  thickness  of  the  transition  layer  is  roughly 
l  Ri  2a  ~  100  lira. 

The  spot  magnetic  diffusivitv  is  -  Vs/vo  —  6.76xl0-2.  For  an  individual 
spot,  the  diffusion  time  is 


ts  =  a2 /£  =  0.25x10-8/6.76x10-2  =  0.037  /xs . 

In  a  distance  X  there  are  XN^2  spots,  but  it  is  wrong  to  assume  that  the 
diffusion  time  for  an  armature  of  side  X  is  simply  the  sum  XN}/2ts,  since  after 
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each  transition  across  a  row  of  spots,  a  weaker  field  is  presented  to  the  next  row. 
The  correct  diffusion  time  is 

A  Y 2 

rx  =  [2aXN}/2)2 /(=—■(.  (7) 

For  the  perfect  contact  surface,  rx  =  rx,  say.  has  the  value  r*x  =  A2//. 
Omitting  the  factor  A/n,  we  get  rx  ~  ctx,  a  physically  obvious  result. 


3.  Pinching  and  heating  at  the  interface 

It  follows  from  the  above,  that  to  replace  the  a-spot  model  by  a  continuum 
representation,  we  need  to  introduce  a  thin  interface  layer,  about  0.1  mm  thick, 
in  which  the  resistivity  is  about  1/e  times  the  resistivity  of  molten  aluminium. 
With  e  =  0.15,  this  amounts  to  a  twenty-fold  increase  over  the  resistivity  of 
the  armature  at  room  temperature.  According  to  the  melt-wave  model,  for  an 
armature  moving  with  a  velocity  V,  the  current  density  has  the  distribution  j  oc 
exp(x/Sr),  where  8r  =  Vr/i^oV)  and  x  is  negative,  being  the  distance  measured 
from  the  melt  wave.  In  the  layer  just  described,  it  would  be  necessary  to  replace 
8t  by  a  number  twenty  times  larger. 

This  raises  the  possibility  that  the  magnetic  field  could  diffuse  rapidly  between 
the  armature  and  rail,  allowing  its  pressure  to  completely  break  the  electrical 
contact  before  the  launch  was  completed.  This  would  result  in  a  pinched  current, 
as  illustrated  in  Fig.  4. 


Figure  4:  Current  pinched  by  magnetic  field 


Consider  the  time  required  for  the  magnetic  field  to  diffuse  half  way  along  an 
armature  of  length  18  mm.  From  (7) 


4x182x10~6 

7TX6.76X10-2 


X0.15  =  910  ^s. 
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But  at  a  speed  of,  say,  V  =  500  ms-1,  the  armature  moves  through  18  mm  in 
36  /is,  i.e.  the  transition  layer  is  renewed  in  this  time.  This  means  that  there  is 
insufficient  time  for  the  diffusion  to  take  place.  However,  if  V  <  20  ms-1,  a  pinch 
will  develop. 


The  ohmic  heating  rate  in  a  spot  is  r]sj 2  per  unit  volume.  Hence  for  the  whole 
surface,  this  rate  is  rjsj2lAs  =  rjsI2l/As.  Per  unit  volume  of  the  corresponding 
continuum,  this  is 


VsP 

ASA 


TlsI2 

eA2 


V 

6 


=  W2  ■ 


In  the  example  given  above,  this  is  twenty  times  greater  than  the  rate  in  solid 
aluminium. 


In  the  steady  state — on  the  continuum  description — the  energy  equation  is 


v-q  =  , 


where  q  is  the  heat  flux  vector.  Integrating  over  a  volume  V,  with  an  element 
dr  =  l  dS,  we  get 

/  V-q  dr  =  [  n-q  dS  =  r/eHj2  IA  , 

JV  JdV 

where  n  is  the  outwards  normal.  Hence  the  jump  in  q  across  the  surface  is 

n  *  N  =  Venf  l , 


or 


Qr  ~  Qa  =  Veffj2  l  ■ 


(8) 


4.  Melt-wave  interpretation  of  the  asperity-plough 
model 


As  just  described,  the  rate  of  melting  in  the  neighbourhood  of  asperities  is  very 
much  greater  than  in  the  ideal  contact  that  formed  the  basis  of  the  previous 
treatment  of  the  melt  wave  model.  By  generating  streaks  of  melted  aluminium 
in  the  transition  layer,  close  to  the  trailing  edge  of  the  armature,  the  c-asperities 
manage  to  cover  a  narrow  region  of  the  armature  with  a  conducting  fluid.  The 
normal  pressure  between  the  rail  and  armature  will  help  to  spread  this  melt 
into  a  continuous  layer.  An  irregular  boundary  will  mark  the  leading  edge  of 
the  collection  of  streaks,  but  shortly  behind  this,  the  streaks  will  run  together, 
forming  a  continuous  layer. 

The  average  position  of  the  boundary  is  a  good  approximation  to  a  melt-wave 
front,  the  depth  h  of  which  is  determined  more  by  the  thickness  of  the  contact 
layer,  than  the  penetration  depth  of  the  current  distribution.  Thus  we  arrive  at  a 
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Figure  5:  The  asperity-plough  model 


model  that  is  a  mixture  of  a-spots  and  a  continuum,  as  shown  in  Fig.  5.  The  rear 
of  the  narrow  conduction  strip  is  terminated  by  an  air  gap  due  to  the  entrainment 
of  the  melt  on  to  the  rail.  Hence  there  will  be  an  ‘entrainment  wave’,  advancing 
along  the  armature,  close  behind  the  effective  melt  wave. 

We  need  an  estimate  for  the  width,  g ,  of  the  region  over  which  the  aluminium 
remains  liquid.  In  terms  of  the  density  density  p,  the  specific  heat  cp  and  the 
thermal  conductivity  K,  the  thermal  diffusivity  is  y  =  Kjpcv.  The  initial  width  of 
the  streaks  of  melted  aluminium  is  2  a  and  the  time  for  the  heat  to  be  (laterally) 
conducted  through  half  this  distance  is  tk  =  a2/x •  During  this  time  the  c-asperity 
will  move  through  a  distance 

g  =  Va2/£ ,  (9) 

relative  to  the  armature.  For  aluminium,  y  =  1.07x10~4,  hence  with  a  = 
0.5xl(T4,  we  get  g  =  0.023H  mm,  i.e.  at  V  =  500  ms-1,  g  &  12  mm.  This  is 
much  larger  than  the  length  scale  for  the  electric  current,  which  means  that  the 
aluminium  strips  will  be  a  liquid  long  enough  to  merge  into  a  continuous,  con¬ 
ducting  layer.  However  we  do  not  expect  the  melted  aluminium  to  fill  the  gap 
between  the  rail  and  armature  surfaces  over  the  whole  width  g.  All  that  the  model 
requires  is  a  narrow  region  of  uniform  conductivity  in  which  electrical  contact  is 
maintained. 

Our  model  thus  closely  resembles  an  ideal  melt  wave.  The  fact  that  the 
melting  initially  occurs  at  discrete  points  on  the  armature,  merely  increases  the 
heating  rate.  But  we  do  require  the  force  between  the  armature  and  rail  to 
remain  high  enough  for  the  melted  aluminium  to  fill  the  gaps  between  the  asperity 
streaks,  as  shown  in  Fig.  5. 

The  derivation  of  the  melt  velocity,  vm  from  the  energy  equation  remains  as 
given  in  [1],  since  the  ragged  leading  edge  has  no  effect  on  averages  taken  across 
the  width  of  the  armature.  The  only  uncertainty  is  in  the  value  to  adopt  for  the 
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ohmic  heating  term.  We  shall  replace  the  ideal  average  resistivity  pa  appearing 
in  the  formula  for  vm  given  in  [1]  by  rf.  Thus 


7T7?*/2  _  G 

2pQw2h  h  ’ 


where 


Q  =  cp(Tm  —  Tq)  +  Lm  , 


(10) 

(11) 


An  increase  in  the  ohmic  heating  term  is  necessary,  since  the  melting  occurs 
largely  at  the  c-asperities.  With  the  ideal  melt  wave  described  in  [1],  we  adopted 
p*  —  fja  =  7.2xl0-8.  The  changes  described  in  §2  imply  that  we  should  use 
the  much  larger  value  p*  =  peK  =  8.5><10_8/e.  A  difficulty  arises  with  e,  which 
is  largely  unknown  and  varies  across  the  region  from  a  small  value  like  0.15  to 
unity.  A  simple  average  like  e  «  0.5  might  be  the  best  one  can  manage,  short  of 
a  very  detailed  study. 

Another  uncertainty  concerns  the  value  to  take  for  h.  Above  we  suggested 
h  =  /,  where  l  is  the  thickness  of  the  contact  layer.  In  [1]  we  showed  from  mass 
conservation  that  the  thickness  of  the  entrained  layer  is  given  by  8m  =  1.1  h vm/V , 
which  from  (10)  is  fortunately  independent  of  h.  We  may  therefore  use 


Sm  =  I- 


ITT]*  I2 

2  pQw2V 


-1 

VI’ 


(12) 


and  measurements  of  8m,  to  determine  p*.  This  will  provide  a  test  of  the  theory. 


5.  Transverse  striations  observed  on  aluminium  rails 

In  Fig.  6  we  reproduce  part  of  a  photograph  taken  at  IAT,  University  of  Austin 
at  Texas,  showing  the  appearance  of  aluminium  rails  following  the  firing  of  an 
aluminium  armature.  The  feature  of  interest  here  is  the  transverse  gouges  in  the 
rail  surface.  From  the  figure  we  find  that  the  wave  length  of  the  ridges  is  about 
1mm.  Apparently  the  melted  aluminium  at  the  contact  surface  has  solidified 
into  waves  with  their  crests  at  right  angles  to  the  motion  of  the  armature.  It  is 
a  reasonable  conjecture  that  these  waves  are  generated  by  an  instability  in  the 
sheared  fluid  layer  separating  the  rail  and  armature.  Vortex  sheets  are  well-known 
to  be  unstable,  a  classic  example  being  that  observed  in  the  wing-tip  streamers 
generated  by  high-flying  aircraft  in  supercooled,  atmospheric  conditions.  Another 
is  the  generation  of  waves  at  the  free  surface  of  a  liquid,  over  which  a  gas  is  flowing. 

The  subject  has  a  long  history,  probably  starting  with  Lord  Rayleigh.  For  a 
recent  treatment  there  is  P.G.  Saffman’s  monograph  [4],  from  which  the  following 
results  are  taken. 
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rigurc  G:  Transverse  striations  in  aluininiuin  rails 


II  tin1  sheai  layer  has  a  thickness  /?,  and  a  velocity  jump  of  {."  (see  Fig.  7) 

waves  of  length  A  =  2tt/ k  have  amplitudes  proportional  to  exp(fr^).  where  tin 
growth  rate  o  is  given  by 

o=±Uk[{l  -  Lkh){l  -  A-//.)]  -  • 

The  validity  of  the  theory  requires  that  h/X  <C  1,  or  hk  «  2tt. 

Applying  the  theory  to  the  sheared  aluininiuin  in  the  railgun.  we  identify 
h  with  h.  the  thickness  of  the  contact  layer  and  V  with  the  velocity  F  of  tin 
armature.  1  he  waves  most  likely  to  appear  are  those  given  by  the  maxirnnn 
giowth  rate',  which  occurs  at  h/t  %  0.G34  and  has  the  value  n  %  O.lTOF/h  It 
follows  from  this  that,  after  the  passage  of  the  armature,  we  should  expect  t« 
see  ’frozen^  waves  of  wave  length  A  =  27rh/().031  =  O.fJlh.  provided  the  tinn 
scale  ,,,  =  ->.8<Sh/l  lor  tlu-ir  growt.h  is  comparable  with  the  time  it  takes  for  tin 
aimatme  to  be  displaced  through  half  a  wave  length,  viz.  r\  =  -I.flGh/F  hhi'- 
rosonaiKv  condition  is  evidently  met. 

If  we  accept  the  estimate  mentioned  in  |j2.  namely  /  =  0.1mm  and  set  h  equal 
to  /.  we  arrive  at  A  %  1mm.  which  appears  to  be  supported  by  Fig.  (i. 

Incident  !y.  the  same  experiment  at  a  later  rime,  shows  rlear  examples  of  the 


Figure  7:  Vortex  sheet  instability 

parallel  striations  caused  by  aluminium  asperities.  These  are  shown  in  Fig.  8, 
p.  12.  The  striations  are  about  0.5mm  apart  and  their  width  is  perhaps  a  half 
of  this  distance,  values  that  correspond  to  a  =  1.25><10_4m  and  Ls  =  5><10-4m, 
which  is  about  twice  the  estimates  given  in  §2  (following  equation  (6)). 
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Figure  8:  Some  evidence  of  streaks  generated  by  a-asperities 
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